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FOREWORD 


Advanced  composites  are  dispersed  heterogeneous  bodies  whose  structure  and 
properties  depend  not  only  on  die  properties  of  the  matrix  and  filler,  but  also  on 
various  kinds  of  interaction  between  the  components,  on  their  mutual 
arrangement  and  volume  ratio  (the  geometric  factor),  on  the  parameters  of  dieir 
IHOcessing  and  production  (the  technological  factor),  etc.  These  very  important 
factors,  whose  change  causes  alteration  of  the  properties  of  a  composite  material, 
must  be  considered  bodi  when  predicting  the  mechanical  behavior,  and  service 
life  of  already  known  composites  as  well  as  when  developing  new  ones.  The 
involved  nature  of  die  interacting  parameters,  makes  theoretical  treatments  of 
these  systems  rather  complex.  This  has  led  to  the  appearance  of  empirical  and 
*  semienqiirical  "structural  and  mechanical"  models  which  are  not  always 
satisfiu:toiy.  Mathematical  modeling  is  one  way  of  solving  technical  problems 
encountered  in  the  development  and  prediction  of  their  mechanical  behavior 
under  diverse  service  conditions.  For  the  construction  of  various  mathematical 
models  it  is  necessary  to  develop  numerical  and  analytical  methods. 

The  modeling  r^iproach  is  especially  important  for  advanced  polymer 
composite  materials.  When  speaking  of  the  merits  and  shortcomings  of  the 
theoretical  approaches  (including  the  molecular  approach)  used  to  describe  the 
viscoelastic  and  mechanical  behavior  of  polymer  matrix,  one  has  to  keep  in  mind 
drat  they  deal  widi  macromolecules  far  from  the  surface.  In  a  composite  (a 
heterogeneous  system),  near  a  filler  surface  one  must  consider  the  interaction  of 
separate  parts  of  a  macromolecule  with  die  surface,  which  could  affect  its 
mobility  and  the  relaxation  time.  The  equation  of  dynamics  of  a  macromolecule  is 
thus  modified  and  becomes  nonlinear.  This  makes  use  of  analytical  methods 
impossible.  Numerical  methods  become  the  best  approach. 

The  theory  of  viscoelasticity  is  based  on  a  function  which  completely 
charactoizes  die  relaxation  spectrum  of  a  material.  In  principle  from  this  one 
pivotal  function  all  other  material  functions  may  be  deduct.  Although  this  seems 
rather  straightforward  in  dieory,  in  jHactice  serious  difficulties  arise,  to  start  with, 
in  obtaining  a  reliable  and  unique  relaxation  spectrum  from  experimental 
measurements. 

To  describe  die  nucromechanical  properties  of  heterogeneous  composites  it  is 
very  important  to  model  their  micromechanical  behavior,  in  particular,  properties 
of  interphase  layers.  For  potymer  composites  the  methods  of  molecular 
mechanics,  molecular  and  Brownian  dynamics  seems  to  be  very  useful.  £>uring 
die  past  few  years,  the  concept  of  a  nanometric  cluster  modeling  of  condensed 
matter,  mainly  of  real  materials  has  been  developed.  However,  when  formulating 
it  quatitatively,  choice  of  the  starting  structure  has  to  be  made. 
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The  objective  of  Ais  workshop  is  to  discuss  the  results  of  recent  studies  and 
experience  accumulated  by  Russian  and  American  scientists  in  the  development 
of  fundamental  theoretical  concepts,  computer  programs,  mediodology  and 
techniques  for  numerical  and  eiqierimental  niodeling  of  structure  and  mechanical 
properties  of  advanced  composite  materials,  in  particular  with  polymer  matrices. 

This  workshop  is  also  expected  to  contribute  to  the  new  contemporary 
concepts  of  "Computer  mechanics  and  Nano-technology". 


US  -  RUSSIAN  WORKSHOP 


ON 

COMPUTER  SYNTHESIS  OF  STRUCTURE  AND 
PROPERTIES  OF  ADVANCED  COMPOSITES 

PROGRAM 

I  Tbuesday,  May  iO,  i^4\ 


9.00-9.30  Registration 


OPENING  REMARKS:  9.30  -  iO.OO  a.m. 
Prof.  Yuri  G.  Yanovsky 
Dr.  Iqbal  Ahmad 
Acad.  Ivan  F.  Obraztsov 


KEYNOTE  LECTURES:  10.00  -  11.30 

CHAIRMAN:  Acad.  Ivan  F. Obraztsov 

10.00  -  10.40  Some  Theoretical  and  Mode!  Approaches  to 

Describing  of  Structure  and  Properties  of  Polymer 
Composites  Materials  -  Yu.G.  YANOVSKY, Institute 
of  Applied  Mechanics,  Russian  Academy  of  Sciences, 
Moscow,  RUSSIA 


10.40  -  11.20  Army  Research  in  Advanced  Polymeric  Composites  - 
A.CROWSON,  U.S.  Army  Research  Office,  Materials 
Science  Division,  USA 


I  11.20  -  12.00  -  Coffee  Break  | 
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Session  A.  TbeoretieaJ  SiamJations 

CHAIRMAN:  Prof.  Yuri  G.Yanovsky 

i2.00  -  i2.40  MatbeautiaU  Method  to  the  Analysis  and  Design  of 

Composite  Materials  -  R.V.KOHN,  New  York  University, USA 

12.40  -  13.10  On  Statistical  Theory  of  Brittle  Crack  Growth  -  A.l. 

MALKIN,  Institute  of  Applied  Mechanics,  Russian  Academy 
of  Sciences,  Moscow,  RUSSIA 

13.10  -  13.40  lU-Posed  Problems  of  Viscoelastic  Media  Equations  from 
Matbemaieal Point  of  View  -  Yu.A.BASISTOV,  Institute  of 
Applied  Mechanics,  Russian  Academy  of  Sciences, Moscow, RUSSIA 


15.00  Lunch 


Session  B.  Nano-Techologicai  Approach 

in  Modeling  of  Composite  Structures 

CHAIRMAN:  Dr.Vladimir  S.Yushchenko 

15.00  -  15.30  Computational  investigation  of  the  Influence  of  the 

Environment  on  Mechanical  Properties  of  Solids.  The 
Cluster  Appmxim^on  -1  .h.GOLMBWih.* ,y  D. 

KH AVRYUCHENKO** ,  A.l . M ALKIN* ,E. A.NIKITINA’and 
V.S. YUSHCHENKO*,  Institute  of  Applied  Mechanics, 

Russian  Academy  of  Sciences,  Moscow,  RUSSIA*, Institute  of 
of  Surface  Chemistry,  Ukrainian  Academy  of  Sciences,  Kiev** 

15.30  -  16.00  Theory  of  Viscoelasticity  and  its  Application  for  Calculation 

of  Relaxation  Properti^  of  Polymer  Composites  •  V.N. 
V.N.POKROVSKY*,  G.V.PYSHNOGRAY**and  Yu.G. 
YANOVSKY***,Economics  and  Statistics  Institute*, Moscow, 

Altai  State  Technical  University,Bamaul**,  Institute  of 
Applied  Mechanics, Russian  Academy  of  Sciences, Moscow, RUSSIA 

16.00  -  16.30  Immibriioa  Modeling  Methods  of  laterpbase  Layers  Properties 

of  Polymer  Composite  Materials  -  S.N.LEORA*,  T.V.IVANOVA**, 
Yu.N.KARNET**,State  University ,St.Petersburg*,Institute  of 
Applied  Mechanics,  Russian  Academy  of  Sciences,Moscow,RUSSlA 
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iS.30  -  i7,00  CriticMl  Periodics  ia  the  Nucleas  Properties  sad 

Problems  of  Composites  Csleulstions  (infoimation  notice) 
V. I. KUZMIN,  State  Institute  of  Radio  Engineering  and 
Automation,  Moscow,  RUSSIA 


DINNER  19.00- 

GEORGIAN  REmAURANT  ’ARAGVI’ 


Wedaesdsy,  Msy  H.  1994 


Session  C.  Computer  Modeling  of  Strength  Stress 
Properties  of  Composites 

CHAIRMAN;  Dr.  A.CROWSON 

9.30  -  10.10  Computer  Synthesis  of  Fiber  Architecture  and  Properties 
of  Textile  Structural  Composites  -  Tsu-Wei  CHOU, 
University  of  Delaware,  USA 

10.10  -  10.40  Asymptotic  Methods  for  Investigation  of  Impact-  and 
Explosion-  Induced  Wave  Processes  -  A.  1.  MALKIN, 
Institute  of  Applied  Mwhanics, Academy  of  Sciences, 

Moscow,  RUSSIA 

10.40  -11.20  Experimentation  and  Modeling  of  Composite  Materials  - 

Shun-chin  CHOU,  U.S.Army  Research  Laboratory, 
Watertown,  USA 

H’TO  -  11.40  Some  Approaches  to  Modeling  of  Composites 

Mechanical  Properties  with  Employment  of  Graphs  Theory 
-  V.U.NOVIKOV,  State  Open  University,  Moscow,  RUSSIA 


11.40  -  12.00  Coffee  Break 


CHAIRMAN:  Dr.  Y.D.S.RAJAPAKSE 

12.00  -  12.30  On  the  Entropy  Damage  Aeeumulation  Model  of  Composite 
Materials  -  S.A.LOURIE,  State  Aviation  Technical 
University,  Moscow,  Russia 


12.30  -  13.10  Coacumnt  Eagineeriag  -  H.T.HAHN,  University 
of  Calofomia,  Los  Angeles,  USA 

13.10  -  13.40  Progress  in  €kunputatioaaI  Model  by  Finite  Element  Method 
of  Viscoelastic  Heterogeneous  Continue  Properties  - 
A.N.  VLASOV*.  V.N.POTAPOV*,A.V.VORONIN**and 
Yu. G.YANOVSKY**, State  Civil  University*and  Institute  of 
Applied  Mechanics**,  Moscow,  RUSSIA 

I  13.40  -  15.00  Lunch  \ 

Session  D.  Advanced  Composites  Properties  and 
Mechanical  Behavior 

CHAIRMAN:  Prof.  Tsu-Wei  CHOU 

13.00  •  13.40  Recent  Advances  in  Mechanics  of  Composites  -  Y.D.S. 

RAJAPAKSE,  Office  of  Naval  Research, Arlington,  USA 

13.40  -  16.20  A  New  Constitutive  Equation  for  Shape  Memory  Composites 
A.A.MOWCHAN,  State  Aviation  University, Moscow, RUSSIA 


16.20  •  17.00  Elastically  Tailored  Composites  Analysis  and  Testing  -  E.A. 

ARMANIOS,  Georgia  Institute  of  Technology,  Georgia,  USA 


17.00  -  18.00  GENERAL  DISCUSSION  &  CONCLUSION 


19.00  -  22.00  CULTURAL  PROGRAM 


8.d0- 

9.30- 

i2.00 

14.30 


9.00- 

10.00 

12.30 

14.00 


1  Thursday,  May,  12,  1994  \ 

19.00  TECHNICAL  VISIT 

Departure:  8.00  from  Hotel  Rossiya 

12.00  Technical  visit  to  Scientific  Industrial  Company 

“Composite’’,  Khat'kovo  (80  kms  outside  of  Moscow) 

■  13.00  Lunch 

■  16.30  Cultural  Program  -  Visit  to  Zagorsk  (Sergiev 

Posad).  Zagorsk  is  the  oldest  scat  of  Russian 
Orthodox  Patriarchy  dating  back  to  11th 
Century. 

Return:  19.00  to  Hotel  Rossiya 


\  Friday,  May  13.  1994  \ 


18.00  TECHNICAL  VISIT 

Departure:  9.00  from  Hotel  Rossiya 


12.30  Technical  visit  to  Moscow  Stala  Aviation 
Technical  University 


13.30  Lunch 

17.00  Technical  visit  to  Scientific-Technical  Centre 

and  Exibition  of  Machine  Engineering  Research 
Institute  by  A.A.BlagotuavoY  of  Russian 
Academy  o  f  Sciences 


Return:  18.00  to  Hotel  Rossiva 
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SOME  THEORETICAL  AND  MODEL  APmOACHES  TO  DESCRIBING 
OF  STRUCTURE  AND  PROPERTIES  OF  POLYMER  COMPOSITES 

MATERIALS 

Yu.G.Yanovsky 

Institute  of  Applied  Mechanics  of  Russian  Academy  of  Sciences 

The  use  of  oi^ganic  (polymer)  matrix  composites  for  different  industrial  application 
actively  started  about  thirty  years  ago.  Comparing  with  traditional  industrial 
materials  there  are  various  advanri^es  that  make  such  composites  more  attractive 
for  use  in  different  kind  of  branches,  in  particular  for  aerospace,  marine, 
automobile  etc.  applications.  Polymer  composites  have  excellent  corrosion 
resistance  which  r^uces  maintenance,  low  thermal  conductivity,  diey  are  non¬ 
magnetic.  Composite  materials  can  be  modified  to  become  radar  absorbing  thus 
enhancing  their  stealth  capability.  Advanced  composites,  featuring  continuous 
reinforcing  fibers  in  polymer  matrixes,  have  demonstrated  their  ability  to  provide 
reliable  weight  saving  structures  for  various  applications.  These  can  offer  weight 
reduction  of  up  to  23-50%  in  comparison  to  aluminum  and  steel. 

It  is  no  longer  a  question  of  whedier  advanced  composites  will  be  used  in  future, 
but  rather  how  to  ^ply  them  more  effectively.  For  modem  polymer  composite 
materials,  which  has  a  history  more  than  25  years,  maximal  experience  has  been 
accumulated  in  aerospace  applications.  Since  aluminum  has  been  used  as  the  base¬ 
line  material  intuition  and  design  judgment  have  been  strongly  influenced  by  the 
aluminum  experience. 

Fundamental  difference  between  metallic,  particularly  aliuninum,  and  polymer 
composite  materials  arise  from:  homogeneity  versus  heterogeneity  and  iimer 
microstructure,  isotropy  versus  macro-  and  microscopic  anisotropy;  elasticity 
versus  pseudo-elasticity  and  microviscoelasticity;  linear  behavior  versus  non¬ 
linear  behavior  etc. 

This  paper  will  focus  on  new  developments  in  the  area  of  computer  synthesis 
(design)  of  reinforce  polymer  composite  materials. 

For  composite  materials  physico-matfaematical  models  provide  basic  information 
to  guide  design. 

The  first  estimates  of  the  elastic  properties  of  composite  materials  fiom  model 
viewpoint  were  undertaken  on  die  basis  of  the  well  known  phenomenological 
models  advanced  by  Voigt  and  Reuss.  These  models  used  the  most  famous  rule  of 
conqxisite  materials  is  the  rule  of  mixtures.  It  state  that  a  composite  property  as  a 
whole  is  a  linear  function  (the  additivity  principle)  of  die  volume  firactions  of  the 
constituents.  In  a  very  simple  variant,  the  behavior  of  a  two-component  composite 
is  modeled  by  using  two  Voigt  or  Reuss  elements  connected  in  pandlel  or  in 
ctmsecutive.  Models  predicting  not  only  elastic  but  also  viscoelastic  properties 
were  the  nmct  stq>  in  this  directitm.  Here  we  must  note  first  of  all  &e  model 
advanced  by  Tak^raiugi. 

Anatysis  of  jdtentMnenological  models  of  the  above  types  and  similar  ones  and 
verification  of  their  adequacy  has  shown  that  at  present  it  is  difficult  to 


recommend  universal  models  for  calculating  die  properties  of  the  composites 
widiin  a  broad  range  of  die  change  of  die  filler  concentration  and  die  conditions  of 
deframation. 

To  be  mtse  predictive  the  {dientHnenological  ^iproaches  and  models  need  to  be 
modified.  A  useful  micromedianical  fiamework  can  be  built  to  model.  If  tbere  is 
no  micnMnechanical  firamework  it  is  nearly  impossible  to  establish  quantitatively 
die  contrilnition  of  the  constituents  and  die  interface.  Empirical  approadi  alone  is 
not  adequate.  A  quantitative  desaipdon  is  important  in  bodi  absolute  and  relative 
terms.  It  is  not  enough  to  say  that  an  interface  is  good  or  bed.  It  is  more  important 
to  know  haw  many  Pascal's  represent  a  good  interfacial  strengdi  and  undmtand 
vrfut  is  a  filler  or  matrix  dominated  behavior.  With  fiameworic  of  rhicromechanics 
one  may  be  rdile  to  clarify  observation  with  quantitative  precision. 

It  is  widely  known  that  the  lack  of  a  stricdy  physical  substantiation  of  the 
phenomenological  models  was  one  of  the  main  causes  of  the  failure  to  describe  the 
properties  of  polymer  composites  observed  experimentally.  From  this  standpoint, 
a  diree-element  structural  model  can  be  considered  more  systematic  and  improved 
in  diis  for  composites.  It  considers  not  only  the  properties  of  the  filler  and  matrix, 
but  also  the  intermediate  layer  spearing  on  the  boundary  of  die  filler  particles 
(interphase  layer)  and  having  specific  properties  that  differ  from  those  of  die 
matrix. 

Direct  experimental  results  obtained  by  birefiingence  and  photoelasticity 
techniques  point  to  die  physical  substantiation  of  die  existence  of  an  interphase 
layer. 

Sununarized  the  stated  above  the  framework  of  the  nanotechnological  ideology 
could  be  written  in  an  appropriate  for  computer  simulation  form  (see  Table). 


Table 

Ideology  of  Program 

"Computer  syndiesis  of  Structure  and  Properties  of  Advanced  Polymer 

Composites" 


Modeling  of  surfaces  and  potentials  of 
interaction  of  separate  phases  of 
heterogeneous  conqiosite  and  (fynamics 
of  its  variation  by  stresses  and 
temperature  (quantum  ^proach) 


Atomic  and  Molecular 
(Nano)-Level 


Macromolecular  Level 


Modeling  of  structure,  mechanical  and 
relaxation  properties  of  interphase 
l^ers  of  conqiosites  by  the  statistic  and 
tmolecular  jdiysics  and  molecular 
medianics  methods. 
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Miax>-Level 


Modeling  of  structure  and  mechanical 
inroperties  of  periodicity  element 


Macro-Level 


Modeling  of  mechanical  properties  of 
structure  representative  element, 
modeling  of  effective  pnqierties  of 
composite  media 


Theoretical  and 
Experimental 
Verification  of 
Modeling  and 
Calculation  Adequacy 


Test  experimental  investigation  and 
mathematical  correction  on  die  base  of 
die  adaptive  program  packet 


The  viscoelastic  behavior  of  a  deformed  polymer  system  is  a  w.^ly  known 
experimental  fact  From  a  physical  viewpoint,  such  behavior  is  generally 
eiqilamed  eidier  by  (he  existence  of  an  intern^  structure,  or  by  the  heterogeneous 
structure  of  the  material,  or  by  both  simultaneously.  The  problem  of  describing  die 
prcqierties  of  such  systems  within  a  broad  range  of  variation  of  die  magnitudes  and 
laws  of  deformation,  the  rates  (or  fi^quencies)  of  deformation  and  other 
parameters  generally  consist  in  determining  die  time  dependences  of  the 
mechanical  diaracteristics,  i.e.  the  material  fimctions  reflecting  die  viscoelastic 
mechanism  of  die  behavior  of  these  systems.  The  ability  of  a  viscoelastic  body  to 
resist  a  mechanical  load  and  react  hereto  adequately  is  determined  by  a  very 
inqrartant  characteristic  of  a  system,  viz.  die  relaxation  time.  For  real  viscoelastic 
media  dut  generally  have  a  ndher  complicated  internal  structure  is  quite  natural  to 
introduce  a  distribution  fimction  of  the  relaxation  time,  viz.  the  relaxation 
spectrum. 

The  material  fimctions  can  be  obtain  firom  eiqierimental  data  rather  simply  but  a 
calculation  of  die  relaxation  spectrum  fiom  above  fimctions  is  inverse  and  ill- 
posed  problem. 

Smoodung  techniques  based  on  regularization  and  the  maximum  entropy  mediods 
have  been  propos^  recently  to  determine  numerically  the  relaxation  spectra  of 
viscoelastic  media  fium  eiqierimentally  obtained  material  fimctions  dirough  the 
solution  of  Fredholm  integral  equaticms  of  die  first  kind.  However,  these  methods 
show  serious  shorteomings  in  q^lications.  We  propose  a  new  and  conqiletely 
different  Hunwteal  tedmiqne  ba^  tm  the  minimax  method. 

A  btge  miodia'  of  fdienomenological  relations  were  proposed  to  describe  the 
viscoeiutic  behavior  of  polymer  bodies.  However  a  phenomenological  desoiption 
does  not  disclose  the  microscopic  nature  of  materials.  The  conadoable  progress 
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•diieved  in  recent  years  in  tbe  field  of  the  fiieotelical  jdiysics  of  polymen  also 
detemined  tbe  achieveme^  in  tbe  structnral  and  mokcuto  dwory  of  pdymers. 
An  c^ialifln  of  die  dynamics  of  a  maoomolecule  in  the  foundation  on  vdudi  a 
mnctnial  dieoey  oi  equilfiviinn  and  noneqnilibrium  lelaxmkm  idmomena  can  be 
contracted  m  a  sin^  molecule  q^ximatum  fra  polymer  systems.  The 
pos^iifify  rtf' imrodiiction  the  single  molecular  qiproadi  is  ba^  in  essence  on  an 
inyntant  experimental  fimt  The  quantities  showing  the  behavira  of  a  polymer 
system  depend  in  a  clear  and  unmnbiguoitt  way  on  die  length  of  a  maoronKdecule. 
When  qicsdcing  of  the  mesits  and  shratcranings  of  die  dieraetical  approadies  used 
to  desraqdioo  die  viscoelastic  bdiavira  of  polymer  we  must  have  in  view  diat  we 
dealt  wiA  macromolecules  fiv  fi'mn  die  surface.  Near  a  surface  one  must  craisider 
the  interaction  of  separate  parts  of  a  macromolecule  widi  the  surfiu».  Here  we 
should  expect  changes  in  die  mobility  of  a  maomnolecule  and  the  relaxation  time. 
The  equation  of  dyruunica  of  a  macrranoiecole  is  modified  and  becomes  non¬ 
linear.  Analytical  appraisal  of  die  results  becomes  inqiossible.  Numerical  method 
of  studying  come  to  die  forefixmt  For  this  purpose,  the  equation  of  dynamics  of  a 
maoomolecule  must  be  written  in  a  form,  fir^,  convenient  for  die  modificatirai 
associated  widi  the  presence  of  surfree  faces,  and,  secondly,  craivenient  for 
pafiMmiag  miinerical  calculations,  because  in  tUs  case  one  hm  to  work  with  a 
aaa  Mncar  system  oi  etpulioos.  One  cm  dnis  state  diat  the  use  of  die  method  of 
Brownian  and  Molecular  (fynamks  in  an  essential  and  impratant  stage  in  studying 
the  merhanism  of  dk  viacoelasddty  of  pedymer  conqwsites  in  boundary  layers. 
Qaantam  chemical  appsoach  to  modd  of  dw  dbesnical  activity,  properties  and 
mterectioB  between  dm  surfiKc  of  a  SUec  and  a  conqiosite  matrix  is  very  inqxxtant 
stage  akag  this  padk  too. 

Qomtmo  chem^  aiproach  has  some  very  impratant  advantages.  Firsdv.  die 
cempatation  of  some  pfaysico-mechanical  chaia^iistics  rai  atom  level,  namely, 
the  ^culated  force  of  atom-atom  interaction,  the  bonds  strengdi  allows  to  obtain 
the  values  of  material  parameters  for  continual  model  and  to  understand  dieir 
sense  and  nature.  Secraidtv.  the  energy,  the  electron  and  space  structure  of 
substance,  srane  intetpfaase  layers  sn^ularity  one  can  calculated  mrae  precisely 
dum  in  classic  molemilar  dyiiamic  version.  As  one  known  molecular  dynamics 
method  are  usually  limited  by  the  quality  at  the  interatomic  pair  potential 
fimetions.  Thinliv.  die  quantum  chemical  computation  allows  to  determine  the 
fiequencies  of  phonon  spectrum  and  the  macroscopic  elastic  constants.  Using  diat 
medrod  tibe  dynamic  model  of  real  substance  can  be  calculated  together  with  its 
GompiesadNlity  and  moduli  of  dasticity. 

It  is  convement  to  enqdoy  die  method  of  finite  elements  to  calculate  die  stressed 
stale  of  die  element  periodicity  and  representative  element  of  media  widi 
periodical  or  non-periodical  structure.  An  object-oriented  finite  elonent 
pcogmmning  seems  to  be  more  praqiecdve  along  diis  path. 

It  has  been  obvious  for  numy  yess  diat  improving  finite  element  code  modnlatity 
demands  a  proper  oiganizmioo  of  die  data.  Maity  audios  have  doie  dieir  bmt  to 
improve  data  management  in  a  Fortran  context  The  adequacy  of  die  concqils  of 
oigectKitientied  progtamming  for  an  efiBcioit  and  nndeidandable  organizatiott  of 
data  in  numerical  analysis  codes  has  only  been  studied  recently.  The 


ialdlllKBtad  principles  at  die  qifdicedoB  of  otgect-oriented  tedinkpirs  to  the 
fliAfc  clement  mediod  were  developed  by  T.Ziniineniieii  (1992).  The  careful 
de^D  of  this  program  resulted  in  a  h^ily  modular,  easily  undmtandaUe  and 
extendable  code.  These  qualities  were  based  on  iticaptitjnii .  aa  olyect  is  a 
imMie  supplied  widi  a  set  a  protected  data  (its  adrflmtes)  tet  allow  it  to  perfimn 
efdiudata  Ots  methodi)  and  to  maaltge  Uteat  operotkms  retfy  antonanondy; 
chaa  iaheritMicg  -  tfci»  varioos  types  Of  objects  (dm  dasses)  are  organized  m  a 
daiple  and  tding  treedike  structare  (the  hierarcby)  vriiich  avoids  code 
dudicatiOBs;  ntm-amk^isriott  -  die  mediods  of  an  obgect  can  be  perfimned  at  any 
piiee  la  die  program,  without  die  coostraming  sequence  of  operations  of  dassicd 
codes. 

An  otgect-miented  apfHoach  was  enqikyed  in  die  new  finite-eleoients  program 
pocket  developed  fijr  cakuladon  of  strengdi-sbress  dasdc,  elasto-plastk  and 
viscoelastic  bdiavior  of  polymer  composites  widi  inter&ce  layns. 

CONCLUSION 

1.  The  structDre-medianical  model  of  heterogeneous  pdymer  conqiorite  widi 
mtetphase  layers  was  enqiloyed. 

2.  Based  on  dus  model  ^  new  qqHoadies  and  computational  methods  of 
atniyhig  of  strengdi<stiess  and  medumical  bdiavior  of  composites  including  the 
nonlinear  viscodasdc  leqxnse  under  deformation  has  been  developed. 

3.  Program  of  calculatirm  crmtains  dm  st^es  vrindi  are  modeling  advanced 
composites  properties  at  die  difGarent  structure  organization  levels:  nano-  and 
mokcular-levd  (quantum  mechanical  metibods);  maoromolecnlar-levd  (Brownian 
and  Mdecular  dynamics  mediods);  micro-level  (object-ruiented  finite  dements 
mediod,  inclodiiig  the  viscoelastic  response),  macro-level  (specialized  pockets  of 
adaptive  computer  programs). 

4.  Present  qifxoadi  is  a  base  of  nano-tedmdogical  princ^e  of  the  heterogeneous 
advanced  polymo’ composites  structme  and  properties  constructing. 
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The  replaceinent  of  heaviCT  metallic  materials  with  advanced  polymeric  matrix 
conposites  of  equal  or  inqnoved  structural  integrity  and  crmsiderable  weight 
savings  offers  significant  (qiportunities  for  the  Army  in  a  variety  of  applications. 
However,  dieir  use  will  require  n  improved  understanding  of  the  relatimiships 
between  processing,  microstructure,  {vqierties  and  performance  derived 
dierefrom.  In  osda  to  achieve  dus  knowledge,  fimdamental  research  is  being 
siqipoited  widiin  die  Materials  Science  and  Engineering  Science  Divisions  at  the 
Army  Research  Office  (ARO).  ^gram  areas  currently  being  addressed  include 
mechanical  bdiavior,  degradation  and  reactivity,  synthesis  and  processing,  smart 
structures,  and  manufacturing  scimice.  In  this  presentation  each  of  diese  areas  will 
be  described  with  examples  of  current  pn^rams  being  siqipoited  at  ARO. 
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The  macroscopic  properties  of  a  contyosite  depend  in  a  subtle  way  upon  the 
microscrqnc  geometry.  It  is  tiierefore  natural  to  lode  for  microstructures  with 
advarrtageous  -  perlups  even  extranal  -  effective  bdiavior.  There  has  been  a  lot  of 
recent  progress  in  this  area,  including  an  identification  of  the  "most  rigid"  and 
"most  crmtyliant"  contyosites  achievable  by  mixing  two  materials  in  specified 
proportion.  Key  tools  include  tiie  "translation  metiiod"  for  bounding  effective 
moduli,  and  "sequential  Umination"  as  a  sdieme  for  generating  optimal 
microstructures. 

A  similar  issue  arises  in  the  modeling  of  shape-memory  alloys.  These  materials 
owe  dieir  special  properties  to  a  martensitic  phase  transfomoation.  Martensite  is 
actually  a  mixture  of  differait  crystallographic  variants,  so  it  may  be  viewed  as  a 
composite  material.  The  microstructure  varies  witii  the  loads  and  boundary 
conations,  driven  by  elastic  energy  minhnization.  Mathenutical  methods  used  in 
die  anatysis  of  extranal  composites  are  also  helping  eiqilain  die  details  of  sh^ie- 
memory  bebavim. 
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ON  STATISTICAL  THEORY  OF  BRITTLE  CRACK  GROWTH 
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INTRODUCTION 

At  tlie  present  tune  there  are  no  commonijr  accq>ted  notions  abont  mechanianu  of 
brittle  and  quaei-brittle  crack  growth  in  eoUde.  It  seems  to  be  associated  not  only 
with  a  host  of  physical  mechanisms  of  growth  but  also  with  an  imperfection  of  the 
existing  system  of  mathematical  modds.  Wthin  the  limits  of  this  report  it  is  im* 
possiUe  to  present  the  grosrth  modds  review  of  any  completeness.  Mention  should 
be  made  of  physically  meaningful  contiaual  modds  (  Barenblatt  et  al),  direct  com¬ 
puter  dmulatinn  by  molecular  dynamics  method  (Yurchenko  and  Shchukin,  Pasldn, 
Dienes  and  many  other  investigators)  and  a  host  of  dislocation  modds  reviewed,  for 
examine,  in  the  monogn4;>h  by  Vladimirov. 

The  vast  majority  of  brittle  crack  models  deals  with  one-dimensional  rectilinear  cracks 
and  treats  the  growth  as  a  sequential  dissociation  of  interatomic  bonds.  Thesemodds 
axe  based  on  a  series  experimentally  found  regularities  such  as  wdl-known  Zhurkov 
longevity  equation.  The  last  is  taken  as  evidence  of  thermofluctuation  nature  of  crack 
growth  at  the  initial  stage  of  fractuse. 

It  has  been  known  that  the  thermofluctuation  growth  of  precritical  cracks  at  brit¬ 
tle  and  quasi-brittle  fracture  shows  the  spasmodic  stochastic  process  characteristics. 
Prolonged  intervals  of  static  behavior  of  crack  ate  alternated  with  fast  moving  of 
the  crack  rip  into  a  new  podrion.  The  distance  between  the  latter  and  the  former 
podrions  of  a  crack  rip  therewith  is  much  more  than  interatomic  distance  in  a  sdid. 
By  this  it  meant  that  notions  of  non-correlated  decohesion  of  interatomic  bonds  are 
not  true. 

ThaoNtial  "ab  initio”  dasaiptioa  of  stoehutie  gtotring  of  brittle  cnek  is  ran  into 
problems  by  the  need  of  taking  into  account  the  cooperative  processes  in  precracked 
aone.  Therefore  the  development  of  phenomenological  theory  based  on  C9q>erimen- 
taUy  observed  regularities  as  wdl  as  on  the  theoretical  reasons  of  extremdy  common 
character,  is  of  particular  interest. 

Here  we  suggest  the  scheme  for  the  construction  of  the  statistical  2-D  theory  of  brittle 
crack  growth.  The  problem  cmuidered  is  concerned  with  several  phenomenolopcal 
modds.  No  considerarion  has  been  given  to  the  detailed  description  of  processes  at 
the  precradced  sane. 

GENERAL  MODEL  FORMULATION 

Let  m  assume  the  crack  at  macrosooinc  scale  can  be  represented  by  the  chain  of  vec¬ 
tors  L  and  hence  specified  by  the  ordered  set  =  {i<i*  =  0,...,n}.  The  dementary 
growth  act  is  the  addition  of  random  vector  to  the  crack  tip  at  the  arbitrary 
random  histaat  of  rime.  In  order  to  formulate  the  kinetic  modd  one  should  construct 
the  probahflity  densi^  of  £»-*  £n+i  trandtion  as  a  fenction  of  process  ^ehistory. 


It  H  dew  that  tiiis  ptobkm  ahonld  be  lolTed  at  tbe  iatennediate  and  micioacoiiic 
•ealea. 

There  are  two  canaea  that  ahonld  be  leeogniaed  tor  the  accidental  nature  of  crack 
(toarth.  That  ia  the  random  arrancement  of  atoma  (or  atructuze  deCecta)  and  the 
aoddental  occupancy  oflocalwed  vibrational  modea  at  tbe  precracked  aooe.  At  given 
ntamk  ocmfignration  the  probaInKty  dennty  of  growth  act  dday  time  may  be  taken 
in  the  cagtonential  Conn.  Thia  aaauTnption  ia  juatified  by  the  rarity  of  growth  eventa 
at  the  time  acale  of  atomic  oadllationa. 

We  ahall  tceept  that  the  detaila  of  the  ptecracked  cone  microatructure  are  not  very 
aignificant,  ao  that  average  delay  time  dependa  mi  a  Umited  number  of  macropa- 
rametera.  It  ia  natural  to  leatrict  our  oonnderatiom  to  the  minimal  aet  of  variaUea 
and  parametera.  Then  ^  average  delay  time  oorreaponding  to  the  defined  growth 
act  r  dependa  only  on  I,  temperature  and  local  atreaaea.  We  auppoae  the  average 
dday  time  to  be  an  Arrenhu  fhnctkin  with  the  activation  energy  linew  -  dependent 
on  atreaaea. 

The  required  invariance  of  the  average  delay  time  in  reapect  to  coordinate  ayatem 
choice  impoaea  the  eaaential  reatrictiona  on  the  form  of  thia  function.  In  the  iaotropic 
medinm  the  average  delay  time  may  be  written  as 

r.(lK)  =  (  ^  -  I  p, 

where  wa  axe  the  components  of  streaatenaor,  I  =  s  and  03(1)  are  aome 

fbnctiana  of  the  adding  microcrack  length,  To  and  Eo  are  acaumed  to  be  constant. 
Expression  (1)  is  a  general  form  of  acalw  r4({|£„)  fulfilling  the  requirements  imposed. 

The  stress  tensor  entering  the  equation  (1)  depends  on  coordinates  of  microcrack 
nndeatkm  point.  In  this  connection  two  widely  used  acenaxioa  of  the  crack  growth 
should  be  recogniied.  In  fine  with  the  first  (microcrack  addition)  scenario  a  micrO' 
crack  arises  firom  some  defoct  at  the  point  I  and  then  is  added  to  the  initial  crack 
tip.  By  the  second  (direct  decoheaion)  scenario  a  microcrack  arises  from  the  initial 
crack  tip  and  then  grows  till  aUqijang  at  the  point  1.  In  that  case  we  — »nmr  the 
microcrack  to  originate  at  the  point  A  =ae  n  with  ae«  I  >.  Here  ae  ia  the  effective 
curvature  at  the  crack  tip. 

Stress  distribution  new  a  crack  tip  can  be  evaluated  under  the  assunqition  that  linew 
fracture  mechanics  is  ^q^calde  to  the  precraAed  aone  description.  For  the  case  of 
isotn^  medium  it  may  be  written 


where  stress  intensity  frctors  (SIF)  dqiend  on  the  initial  crack  configuration.  Ex- 
presaion  (2)  cocteaponda  to  aeenario  of  direct  decdieafam  at  the  crack  tip.  However, 
H  ia  eaigr  to  chedr  that  (with  minw  reservation  }  general,  form  doea  niot 

dq^end  on  the  scenario  of  crack  qrawth.  NcAice  thad  substitntimi  of  ae  doet  not  add 
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the  munber  of  empirical  paninetera  m  m  and  Vj(t)  enter  into  eq.  (1)  in  the  fimn 

It  ia  b^ered  that  poatulated  eq.  (1)  cmreaponda  to  some  ^proodmation  of  the 
average  fidd  type.  If  this  ia  the  caae  the  reasonable  interpretation  of  the  equation  (1) 
may  be  provided.  Thia  interpretation  ia  based  im  the  supposition  that  local  stresses 
at  the  point  of  microcrack  nudeation  have  a  special  form 


+  qi{l)nink<Tik  (3) 

where  om  >•  mactostress  tensor  (2)  (ojs  proposed  to  be  the  result  of  the  averag¬ 
ing  over  the  statistical  ensemble),  fonctions  qi(I)  and  qi(l)  account  for  the  excessive 
stresses  caused  by  the  structural  defects.  These  functions  should  be  defined  at  the 
intermediate  scale.  The  stress-induced  reduction  of  the  activation  energy  of  microc¬ 
rack  nudeation  is  defined  by  activation  vdumes  and  which  does  not  depend 
on  1.  Activation  vdumes  are  proposed  to  be  the  material  constsmts  and  should  be 
calculated  at  microscale.  It  is  easy  to  check  that  each  vj{l)  in  eq.  (1)  is  the  bilinear 
form  of  smd  yk(0‘  Needless  to  say  that  the  latter  functions  are  difficult  if  not 
impossible  to  obtain  theoretically. 

We  shall  introduce  the  state  density  function  0(1).  This  function  represents  the 
number  of  microscopic  distinguishable  states  which  lead  to  the  dementary  growth 
act  for  a  vector  1.  Let  us  consider  the  virtual  set  of  the  growth  acts.  Clearly  the 
growth  act  to  be  realised  has  a  minimal  random  delay  time  t(i).  Therefore  non- 
conditkmal  distribution  of  growth  act  dday  time  is  defined  as  the  distribution  of 
random  value 


t„  =  min  f(i) 

{0  (4) 

where  the  contribution  from  t(i)  should  be  taken  into  account  with  a  weight  0(1). 
Rsmdom  value  (4)  probabifity  distribution  is  also  exponential.  The  average  value 
<  tm  >  i*  determined  from  the  expression 


^  Im  > 


o(i)df 

r4ii\Cn) 


(5) 


where  rtiJ^Cn)  is  given  by  expressions  (1),(2).  As  a  result,  a  posteriori  distribution 
of  random  vectors  1  is 


_  «(Q  ... 

"  ZrsdK) 

The  rdatioBS  presented  oam|detdy  determine  a  set  of  crack  growth  models  when  the 
fimctioos  Vj(l)  and  state  dennty  function  0(1)  are  prescribed. 
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N«^  that  if  ttate  density  fuactioa  does  not  depend  on  stresses,  then  it  wiO  be  the 
fonctioD  on  tnictocrsck  length  but  not  the  direction,  0  =  n(l).  In  that  case  (!({) 
represents  only  the  structure  properties  of  the  solid  at  intermediate  scale.  The  last 
conclunon  is  true  fat  isotro{Hc  media. 

SOME  PROPERTIES  OF  ADDING  MICROCRACKS  DISTRIBUTION 


Geometric  characteristics  of  2-D  cracks  contour  are  defined  by  function  (6).  The 
analysis  of  this  function  leads  to  some  restrictions  on  In  particular,  the  typical 
crack  contour  form  may  by  obtained  with  the  proviso  that 


ei(l)  >  0,  -2ei(I)  <  »,(/)  <  -vi(/)/4  (7) 

Let  us  assume  that  the  initial  crack  is  close  to  rectilinear  mode  I  one,  so  that 
KttjKx  1.  If  inequalities  (7)  are  fulfilled,  then  the  most  probable  direction  of 
growth  will  be  close  to  the  direction  of  the  initial  crack.  The  conditional  probability 
dennty  of  growth  direction  becomes  the  simple  form 


tD(d|l,£n)  = 


y/2ic  <Sfi> 


:e  !<«'> 


<  id*  >S= 


ikTy/ZKK 

V.Kl{Cn) 


where  o.  «  2vi(l)  —  vj(/)  Basing  <m  the  need  to  obtain  the  correct  expression  for 
specimen  lifetime  (i.e.  Zhurkov  kmgevity  equation)  one  can  conclude  that  microcrack 
length  dependence  of  vy  is  bound  to  be  smooth.  Mote  precisely,  vi  and  vj  variation 
with  microcrack  length  must  be  small  for  the  range  contributing  significantly  to  the 
integral  (5).  An  alternate  possibility  consists  of  a  special  definition  of  functions  vj(l) 
and  n(/).  This  posnbility  seems  to  be  highly  improbable. 

By  these  means  the  vj  dependence  on  the  adding  microcrack  length  should  be  ne¬ 
glected  at  least  at  the  initial  stage  of  crack  growth.  It  is  therefore  concluded  that 
the  mictocrack  length  and  direction  are  independent  random  values.  As  a  result,  the 
distribution  density  of  mictocrack  length  (6)  may  by  rewritten  as 


w(l|£»)  =  w(l)  =  I"  (9) 

and  does  not  depend  on  process  prdiistory.  The  first  order  correction  to  the  most 
probable  microcrack  length  I,  is  represented  by  expression 

».  Af*  Ki  ^1 

-  (2ir4)*n»(Mfcr  ®  di  I, =/. 


where  v+  =  2vi  +  oj  and  derivation  dv+/dl  iq>pears  to  be  always  positive.  The 
condition  fm  fomtola  (9)  validity  may  be  defined  as  Sl/le  <  1. 

GENERAL  KINETIC  EQUATION 
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Within  the  Hxnitt  of  ttatuticel  approech  it  ia  natural  to  describe  the  crack  growth 
by  the  probability  denaity  of  transition  from  the  initial  state  Co  to  the  state  iCn  m  a 
time  t .  We  shall  denote  this  function  by  Based  on  the  suppositions  above 

it  is  not  difficult  to  obtain  the  basic  equation  in  V{Cn,t) 


n(i,)(i-<»n) 

T4(r«|£,_x) 


ViCn-ut) 


The  initial  condition  for  eq.  (10)  is 


V(£„,0)  =  SonS(Co,C^^'>) 


n  =  0,l,...  (10) 


(11) 


Here  is  Kronecker  delta,  S{Co,C^^)  is  delta-fimction,  £0°^  is  the  prescribed  initial 
crack  contour. 

lb  solve  the  initial  value  problem  (10),  (11),  SIP  at  every  growth  step  should  be 
known  .  Analytical  approach  is  possible  only  in  the  small-angle  approximation.  Here 
we  shall  restrict  our  consideration  to  the  more  simple  case  of  rectilinear  crack  growth. 

RECTILINEAR  GROWTH  APPROXIMATION 


Let  us  consider  near-the-suriace  crack  in  a  semi-infinite  specimen  (i.e.  in  a  half¬ 
plane).  One-dimensional  approximation  of  rectilinear  growth  corresponds  to  the  Umit 
Kii  — »  0  ,  <  6B^  >—*  0.  Consequently,  the  conditional  distribution  (8)  transforms 
into  delta-functi :>n,  and  crack  configuration  is  characterised  by  the  unique  parameter 
•  the  length  L.  If  the  crack  contour  is  normal  to  the  specimen  boundary  then  SIP 
Jf/ a;  l,12V»i<r«,. 

Distribution  density  of  rectilinear  crack  length  should  be  introduced  as  follows: 

Vr(L,  t)  =  f;  /  dCndL^S(L  -  nyViCn,  t)  (12) 

n=0  ■'  <=0 

where  dCn  =  Hilso  h  =  Co  w  a  vector  assigning  the  initial  crack  tip  location. 
In  order  to  obtain  the  equation  in  l^(L,t),  one  should  multiply  eq.  (10)  by  S{L  — 
integrate  the  product  over  dCndL^  and  then  summarise  the  result  over  all 
n  considering  V’(£„,t)  =  0  for  n  <  0.  In  consequence  of  this  procedure  we  have 

O’) 

-I-  £  dVw(L  -  £')e^v^ 

with  the  initial  condition 

v;(£,0)  =  «+(!- is) 

Here  the  dimensionless  variable  r  and  parameter  $  are  pven  by 
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•ad 


^  l,12v+iroa  [L^ 

■  kT  V2I; 

Nomerical  Kdotkm  at  eq.  (13)  ahoiild  not  present  any  essentia]  difficulties.  As  for  the 
analytical  studies,  there  it  a  more  convenient  way  of  treatment.  This  way  is  based 
on  the  so-called  difEurion  model. 

DIFFUSION  MODEL  OF  RECTILINEAR  GROWTH 


Elquation  (13)  can  be  reduced  approximately  to  the  more  readily  solved  partial  dif¬ 
ferential  equation.  This  equation  may  be  obtained  in  the  usual  way  when  the  major 
contribution  to  the  int^al  term  in  eq.  (13)  is  given  by  a  small  range  near  the  upper 
limit.  It  can  be  shown  that  the  necessary  condition  for  the  diffusion  model  validity 

</><  — - - 


V+ffoa 


The  diffusion  model  equation  of  rectilinear  crack  growth  takes  the  form 


dr 


=  -<!> 


<e»> 

2 


(14) 


-1/4. 


It  is  clear  that  eq.(14)  represents  the  crack  growth  process  at  the  crude  scales  of 
length  and  time.  The  results  of  the  diffusion  model  are  meaningful  when  the  growth 
observation  time  interval  is  considerably  larger  than  the  growth  elementary  act  ex¬ 
pectation. 

GROWTH  KINETICS  IN  THE  DIFFUSION  MODEL 


It  is  convenient  to  conduct  growth  kinetics  consideration  by  the  use  of  the  equation 
in  the  integral  distribution  function  U{Vr  =  dV/dL).  Substitution  x  =  —1  +  £/Lo  > 
r  — »<  I  >  r/Lo  leads  to  the  fcdlowing  form  of  this  equation: 

where  D  =<  P  >  /2<l>  L„. 

The  initial  and  boundary  conditions  for  U  are 


U(a,0)  =  1,  *  >  0 

J7(0,t)  =  0,  U(oo,t)  =  1,  T  >  0  (16) 

The  s<dution  of  boundary-value  problem  (16),  (17)  can  be  obtained  by  numerical 
simulation  mily.  However,  if  D  >  1,  the  approadmate  solution  is  constractible.  This 
stdution  is  represented  by  the  asymptotic  expansion  over  small  D  which  is  uniformly 
suitable  at  a  wide  range  of  independent  variables. 
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Some  notknu  of  the  crack  growth  kiaetica  can  be  gained  from  the  first  atatiatical 
momenta.  The  cakoUtiona  were  carried  out  bjr  the  uae  of  approximate  solutkm.  In 
the  moet  interesting  case  ^  >  1,  for  the  dimenaionlesa  average  crack  length  1+  <  x  > 
and  dispersion  <  >  we  obtained  the  following  expiessiotts 


where 


1+  <  *  > 

<  fx*  > 


-  t) 
2 


4De-»<»  .s/9(r.-T) 


+  ••• 


T, 


2P 


(17) 


Expressions  (17)  are  the  dominant  parts  of  the  moments  expansion  over  smaU  D. 

It  should  be  recalled  that  the  considered  modd  is  suitable  for  initial  stage  of  the  crack 
growth.  Therefore,  expressions  (17)  do  not  describe  the  growth  at  the  r  -*  r,  limit. 
The  crack  tip  tends  to  infinity  at  the  finite  instant  of  time  is  due  to  the  neglecting  of 
growth  elementary  act  duration. 

It  is  of  interest  that  the  crack  length  dispersion  shows  an  extremely  rapid  increasing. 
The  ratio  of  distribution  half-width  to  average  length  diverges  as  (r,  —  r )  ln~^  (r.  —  r ). 

CONCLUSION 


Here  we  have  tried  to  formulate  the  phenomendogical  approach  to  the  statistical 
description  of  brittle  crack  growth.  We  discussed  the  limiting  case  of  an  extremely 
localised  growth.  Furthermore,  our  consideration  was  restricted  to  the  minimal  num¬ 
ber  of  statistical  degrees  of  freedom.  The  growth  elementary  act  was  assumed  to  be 
defined  by  geometrical  parameters  of  adding  microcrack  only. 

Apparently,  the  suppositions  used  are  rather  limiting  for  real  crack  growth  processes. 
However,  the  number  of  generalisations  is  possible  by  the  method  ^ven  above  without 
essential  modifications.  This  is  particularly  true  for  substitution  of  added  statistical 
variables  which  specify  the  growth  elementary  act.  It  is  desirable  that  the  method 
which  is  under  discussion,  can  be  used  for  the  construction  of  more  realistic  statistical 
modek. 


i  ' 
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The  linear  thecny  of  viscoelasticity  is  based  cm  one  function  i^ch  completely 
characterizes  die  matoial,  die  relaxation  spectnim.  In  princ^le,  frmn  this  one 
pivotal  functkm  all  odier  matmial  functions  nu^r  be  deduced.  Althou^  this  seems 
radier  straightforward  in  dieoiy,  in  inactice  serious  difficulties  arise,  to  start  widi 
in  obtaining  a  rdiable  and  unique  relaxation  spectnim  from  experimental 
measnremoits. 

The  difficulties  encountered  are  two  fold.  Firstly,  considerable  noise  is  superposed 
on  die  discrete  values  measured  depending  on  die  mediod,  and  secmidfy  die 
domain  of  die  disoete  values  is  certahily  restricted  in  any  eiqieriment.  Fuidier,  the 
inverse  problem  of  finding  the  relaxation  spectrum,  die  kernel  of  a  Fredholm 
equation  of  die  first  land,  from  measured  data  is  an  ill-posed  problem.  Inverting 
dtis  int^ral  equatitm  is  not  an  easy  task  because  the  solution  to  die  inverse 
problem  is  not  unique.  It  can  be  inade  unique  only  by  introducing  additional 
assunqitions  ot  additional  a  priori  information.  Recendy,  Tikhonov  regularization 
and  entrcqiy  fimctioiud  methods  have  been  qiplied  to  the  inverse  problem  of 
finding  tite  kernel  of  die  Fredholm  equation.  Both  mediods  allow  unique 
determinatioa  of  the  konel  by  introducing  a  priori  error  estimates  and  optimizing  a 
r^jularization  functimial.  fritroducing  error  estimates  and  assuming  a  certain  error 
distribution  makes  the  solution  unique,  but  not  necessarily  the  correct  solution  as  it 
depends  on  die  a  primi  estimates.  Furdier,  in  any  experimental  data,  in  addition  to 
die  inherent  noise,  there  is  mnbedded  another  set  of  eiqierimental  values  vdiich  are 
not  accounted  for  by  any  error  estimate  and  distribution,  called  outliers.  The 
answer  to  die  questitm  as  to  how  large  is  diis  set  depends  on  many  factors.  But  it  is 
certainly  not  negligible  by  any  means.  Lasdy,  we  remark  diat  if  the  data  has  been 
obtained  in  a  restricted  region,  which  it  is  almost  always,  extending  the  relaxation 
spectrum  deduced  from  it  by  i^iatever  to  a  larger  range,  which  is  always  necessary 
in  practice,  is  a  moot  question  by  definition  if  the  method  used  to  obtain  the 
qiectrum  in  the  limited  range  has  shortcomings  to  begin  with. 

The  new  algoridun  based  on  die  minimax  mediod  we  develop  does  not  require  a 
priori  earn  estimates  and  die  set  of  oudiers  are  taken  into  account.  We  construct 
the  algoridun,  give  several  hypothetical  examples  to  test  die  reliability  and 
exactness  of  die  method,  and  finally  end  die  piqier  widi  a  concrete  example  of 
relaxation  qiectrum  we  determine  for  a  polymeric  liquid. 
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INTRODUCTION 

The  phenomenon  of  the  sharp  decrease  of  the  strength  of  soKds  in  some  enrironments 
(specific  for  each  material)  is  well-known  and  finds  many  ^plications  in  various 
branches  of  en^eering.  Siometimes  (especially  for  liquid  metal  embrittlement)  the 
rednctimi  of  strength  may  be  some  tens  or  hundred  times.  This  effect  appears 
immediate^  after  the  contact  of  stdid  with  such  active  media  and  is  quite  different 
from  usual  corrosion  or  stress-corrosion.  Obviously  the  active  medium  participates  in 
the  process  of  the  fracture  -  if  the  acUve  medium  is  removed,  the  strength  returns  to 
its  initial  value.  This  phenomenon  has  been  studied  for  many  years  from  pioneering 
work  of  P.A.Rehlnnder  (1028)  who  related  the  effect  with  the  reduction  of  snrfrce 
free  energy  of  sdid  in  contact  with  the  active  medium  and  the  decrease  in  the  work 
of  fracture  because  the  fracture  always  is  the  creation  of  new  surfaces.  This  idea 
has  found  various  confirmations  at  least  as  a  necessary  condition  of  the  phenomenon 
that  is  often  named  as  the  Rehbinder  effect  (RE).  Such  thermodynamic  approach 
may  be  used  as  the  bans  at  the  macrosco^c  analysis  of  the  RE,  but  it  leaves  aside 
the  mechanism  of  the  phenomenon. 

There  are  various  hypothesis  about  the  atom-scale  mechanism  of  the  RE,  but  none  of 
them  is  able  to  include  all  cases  of  its  manifestation.  In  last  decade  the  understanding 
speared  that  the  RE  is  the  complex  phenomenon  that  may  be  caused  by  various 
processes  at  the  boundary  between  the  sdid  and  the  active  medium  and  it  is  more 
comet  to  speak  about  the  possible  specific  mechanisms  of  the  fracture  in  active 
medium  and  thw  rede  in  total  influence  of  environment  on  the  mechanical  behavior 
of  the  s<did. 

It  is  commonly  supposed  that  the  most  ’’direct”  and  obvious  mechanism  RE  is 
decohesion  -  reduction  of  the  strength  of  interatomic  bonds  in  the  surface  layer  of 
the  solid  interacting  with  active  me<fium  (and  active  media  are  those  that  weaken 
interactions  in  s<dids).  The  quantum-chemical  calculations  of  various  level  seem  to 
be  the  only  way  of  theoretical  investigation  of  this  mechanism. 

The  omnplex  nature  ai  the  effect  makes  it  very  difficult  to  use  first-i»inciples 
iqtproadies  to  study  the  processes  of  deformation  and  fracture.  There  are  two  ways 
to  overcome  this  problem:  (1)  to  choose  the  qrstem  sim]de  enough  that  really  may 
be  considered  frmn  the  first  principles  and  (2)  to  use  some  of  semiempirical  methods 
that  allow  to  deal  with  larger  systems  and  to  emusder  mote  or  less  large  dusters  of 
atoms. 
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it  ■wmi  tliat  OKMt  if  not  all  of  tlie  kaown  qaantam-diieiDical  calcuUtioni  have  been 
den*  hf  the  &et  method.  In  particnUi:,  methodkai  *Ak  «mt*o^  cakmlatioBa  are 
ftiWBed  (  V.8.  Ynahchcoko,  BJ).  Shdinkia  et  aL)  for  rather  aman  syetema  with  the 
C  —  C  and  C  »C  booda. 

TUa  wagr  ia  aomewhat  anbatantiated  at  leaat  for  quaEtative  or,  at  beat,  aemi- 
qoaatitative  deacriptiaa  of  the  efleet  ia  eoEda  with  coralent  (or  ^edondaaatlj 
oofaleat)  boada.  Ia  that  caae  the  piopertiea  of  the  iateractiona  are  detenniaed  bjr  the 
iateraetiag  atoma  aad  their  neareet  neighbors.  A  different  atnation  ia  reahaed  ia  the 
metaWr  aaKda  Aaractenaed  by  the  oooqdetelgr  aon-local  interatomic  iateractiona.  It 
ia  apparent  that  metallic  aolida  ahonld  not  be  atndied  with  the  haage  of  amaU-acale 
atomic  ayatema.  Becanae  ci  thia  the  aecood  wi^  ia  preferable  over  the  metallic  and 
ionic  aolida.  However  even  ia  the  caae  <d  covalent  aolida  it  k  queationable  whether  the 
leanlta  obtained  fm  amaQ  molecuka  are  uaefnl  in  the  deaciiption  of  fract  nre  proceaaea 
in  the  crack  tq>.  The  point  k  that  even  a  tini^  act  of  adaorption  may  produce  the 
agnificaat  atomic  rearrangement  in  a  moderatdy  large  domain  cA  the  atreaaed  solid. 
Thk  kind  oi  inatabilitiea  aeema  to  be  impmtaat  in  fracture  proceaaea. 

By  these  means  the  aemiempirical  cluster  methods  are  required  few  quantum- Aemical 
investigationa  RE  despite  their  obvious  pitfaUa.  Cluster  ^proach  has  the 
advantage  that  the  adsorption-induced  rearrangement  of  the  precracked  aone  can 
be  calrnlated  directly.  Here  we  aemiempirical  AMI,  MNDO/H  and  NDDO 
methods  aUe  to  operate  with  500-orbital  systems  with  comidete  optimisation  of 
atomic  structure  for  the  investigation  of  the  adsorption-mduced  biealmge  of  stressed 
hydrocarbon  chains,  dusters  of  aluminium  of  different  rise  as  wdl  as  with  the  dusters 
on  immovaUe  substrate  and  with  bulk  silica  systems. 

To  test  the  validity  of  these  methods  for  the  proUem  under  consideration  the 
comparison  between  ”ai  instto**  calculation  results,  and  semiemprical  ones  was 
performed  for  the  simplest  systems. 

THE  CALCULATION  METHOD 

The  quantum  Aemical  methods  MNDO/H  and  AMI  applicable  for  a—  aad  p- 
dements  have  been  offered  by  Dewar  in  1985.  These  methods  are  based  on 
sdfeanaktent  aemiempirical  apptoadi  giving  the  predrion  comparaUe  with  the  ”ab 
Msstso*  method  STO  —  4  —  310F  bat  10*  times  quicker.  We  suppoae  these  methods 
(  with  the  modifications  required)  to  be  the  most  convenient  tot  the  RE  studies. 
The  modifications  indude  "energetic”  parametrisathm,  the  opportunity  oi  the  d- 
orlritak  involving  for  transition  dements  and  the  algorythm  tot  the  deformed  states 
simulation. 

Tlie  suggested  duster  approrimatiou  makes  it  possible  to  folfiQ  the  total  optimisatum 
of  the  spaee  structure  of  systems  oontaiaing  up  to  ISO  -  300  atmns  using  PC 
eompntecs  based  on  INTEL  480/487. 

The  computer  dmuiarion  of  RE  in  the  duster  i^noach  includes  two  main  stages. 
The  pn&niaaiy  study  involves  the  moddHng  of  the  systems  with  various  initial 
eoufignsadoM  fo  oeder  to  reveal  the  aet  of  low-energy  duster  structures.  The  latter 
k  very  impestant  for  rimuhtiou  of  the  riow  adiabatic  decohesion  proeeeaes  because 
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il  k  aatiml  to  axpact  Um  fiaito  cnafiguratioM  to  be  the  moct  etable  low-energy  me*. 

One  of  the  ntkn  dHBcohiee  of  the  dnetcr  i^peoech  cooms  from  the  i»oUem  of  the 
heee  doeter  dtoiee.  Wb  ehoold  eeaidi  for  the  •malleBt  miteUe  ccmfignration  in  mder 
to  simidify  the  cefcoletione.  At  the  Mine  time  the  doeter  moet  be  large  enough, 
capable  not  only  to  deactibe  the  propertiee  of  the  eui&ee  but  aleo  to  modd  all  the 
poeaMt  adaocptiaa  dtea  for  difforent  types  of  adeorbatea.  In  acoorduioe  with  the 
theory  of  the  "aetioe  anaemfales”  the  problem  is  to  rereal  how  maiqr  atmns  should  be 
ifthn  comidflnticMia 

‘Hm  problem  of  the  duster  diaiee  differs  for  different  types  of  systems.  For  the 
eosalent  affiea  system  the  duster  dioold  contain  the  combination  of  seeeral  siUccm 
•  oxygen  tetrahedra  surroundiag  the  Si  —  0  —  Si  bond.  For  the  metalHc  systems 
the  duster  can  be  choeen  <m  the  base  of  the  bulk  crystal  structure  or  the  definite 
crystaBographic  face.  In  this  case  the  number  of  atoms  of  the  upper  layer  necessary  to 
be  taken  into  account,  depends  on  the  migin  of  the  metal  and  the  adsorbate  as  wdl. 
It  is  also  not  clestr  how  many  layers  take  part  in  the  adsorption  and  deformation 
process  especially  for  the  metals  tot  which  the  strong  adsorbate  induced  surface 
reconstruction  is  observed. 

The  second  stage  indudes  the  study  of  the  influence  of  adsorbates  on  the  defiormatkm 
properties  of  the  qrstems.  In  <ader  to  solve  the  main  problem,  the  basic  method  was 
adapted  for  calculating  of  different  types  of  medmnkal  deformations  of  dusters  such 
as  uniaxial  defiormatkm  and  bending. 

The  computational  coqieriment  consists  of  sequential  step-by-step  deformation  from 
the  initial  stable  state  to  rupture.  The  elementary  deformation  step  is  choeen 
acoisrding  to  the  value  of  interatomic  distances  in  the  non-deformed  duster.  The 
proper  step  length  is  found  to  be  10~*...10~*  of  the  bond  length.  The  comidete 
optimisation  of  space  structure  is  perfcwmed  at  each  step. 

It  is  dear  that  such  a  method  is  able  to  describe  a  sequence  of  staUe  states  but  not 
the  kinetics  of  rupture.  In  order  to  obtain  the  kinetic  characteristics  one  should  turn 
to  wdl-known  models  of  the  reaction  rate  theory.  As  we  are  able  to  calculate  the 
vibrational  spectra  and  saddle  points  parameters  the  kinetic  characteristics  can  be 
estimated. 

The  method  described  wh  iq>]died  to  the  investigatioa  of  the  mechanical  properties  of 
hydrocarbon  dmins,  bulk  silica  systems  and  aluminium  dusters  under  the  adsmption 
of  several  water  molecules,  proton,  hydroxjdt  l^draxonium  kms. 

THE  APPUCATION  TO  HYDROCARBON  CHAINS 

In  order  to  compare  the  results  obtained  for  RE  by  the  method  under  discussion 
with  those  obtained  by  eh  tntfio  caknlatwins  we  have  reproduced  the  energy  and 
Mtaratomk  focees  calculations  of  double  C  -C  bond  in  ethylene.  The  influence  of 
the  hydrogen  cation  on  the  bond  weakening  was  also  studied.  The  results  obtained 
by  the  semiempiiical  and  eh  iniUo  calcnlations  are  in  good  agreement. 

The  sdenM  of  the  moddflng  of  *  dsfermation  types  (the  niii»Ti«l  dcformatkm 

and  the  bending  )  was  tested  on  a  butane  molecule.  The  dependence  of  the  space 
and  siect  runic  stencture  and  tbs  heat  of  formation  on  mecJtanical  deformation  of 


tike  wwlenile  iMte  obtained  ae  weD  ae  tbe  eSicatioBal  q^ectra  far  the  free 

and  stntched  heptane  mobcnk.  The  «d>taiaed  loir  freqiaency  shift  far  the  stretched 
molecnle  fa  in  fnaitaftee  aeeerdaaee  with  e^ecfanental  data  far  po^rmets. 

THE  APPUCATION  TO  SIUCA  SYSTEMS 

One  farther  exsantde  of  the  application  <rf  the  method  under  discnssion  to  the  coralent 
spetems  fa  the  study  the  morhanical  properties  at  siBca  dfaxide  in  active  media. 
The  mediaiiical  stretdi  of  the  5*  -  O  —  Si  bond  was  examined  in  vaeunm  and  in  the 
preaenoe  of  several  reagents  (water,  the  hydraranium  ion,  some  snrface-active  organic 
molecniee,  etc.).  The  clnater  containing  two  aKca-oxygen  tetrahedra  «no«i«Jltng  the 
active  center  of  aiKca  dioxide  was  chosen.  The  alteration  of  the  heats  of  farmation, 
the  changes  in  the  stractnral  and  decttmuc  configurations  of  the  clusters,  and  the 
farce  characteristics  of  deformed  systems  were  obtained  (Fig.l).  The  hydroxonium 
ion  was  shown  to  have  the  greatest  inflnence  on  the  rupture  of  the  Si  —  O  —  Si  b<md. 
It  may  be  conduded  that  the  bond  decoheaon  force  in  the  slightly  acidic  medium 
fa  reduced  ai^eximately  three  times  and  the  work  of  rupture  -  almost  ten  times 
being  compared  with  that  for  the  free  cluster.  The  dose  results  were  obtained  in  the 
presence  of  the  surface  active  substances  (SAS). 
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THE  APPUCATION  TO  IfETALUC  SYSTEMS 

Hia  ptdimiiury  ttag*  Hie  iimatigetiaa  iaclndiag  the  adectioa  of  the  defter  nae 
ead  ftructun  beoanea  for  nMtalHc  tTstema  very  important.  The  duater  should  be 
large  enough  to  poaaeff  the  prr^erties  ( ioniiation  potential,  the  lace  crystaDogn^hy, 
etc.)  doae  to  the  bulk  metal  onea.  On  the  other  hand  the  system  should  be  sufficiently 
small  to  be  accessible  tot  serial  calculations. 

Ibtal  optmuaatkm  ot  dnaters  containing  from  2  to  26  di  atoms  give  us  an  opportunity 
to  make  certain  condnaiona: 

(1)  several  isomers  can  be  obtained  for  almost  all  of  the  dusters  of  each  nae; 

(2)  the  most  low-energy  modifications  of  small  aluminium  dusters,  containing  from 
2  to  7  atoms,  are  planner; 

(3)  the  stable  dusters  containing  from  8  to  11  atoms  are  transitional  between  planar 
and  3-dimensioaal  stmctnies; 

(4)  the  bulk  structures  are  obtained  Im  the  dusters  with  the  number  of  atoms  more 
than  12; 

(5)  the  duster  of  AJi$  possess  the  lowest  heat  of  formation  and  the  total  energy 
rdated  to  one  atom  among  all  the  dusters  from  Alj  to  Aljt,  it  appears  to  be 
a  stable  system  of  tetrahedra  connected  in  twos  by  sides  and  by  two  common 
tops,  the  lqrbridisati<m  the  orbitals  of  the  central  atom  u  sp*; 

(6)  all  the  bulk  dusters  seem  to  be  built  from  stable  Ali$  block  by  fiHing  the 
successive  coordination  sheila  by  Al  atoms  forming  the  additional  tedrahedra. 

Accordiiig  to  the  reasons  above  we  have  chosen  the  duster  as  the  main  object  for 
investigation  of  RE  for  aluminium  systems.  AI21  looks  like  two  blocks  having  a 
ctmunon  quadrangular  nde. 

Under  the  influence  of  external  stresses  the  Aljs-block  is  deformed,  the  reconstruction 
of  all  the  bonds  occurs.  We  have  found  several  sones  of  elastidty  separated  by 
the  "plastic”  sones  <d  two  types:  the  sharp  reduction  of  stresses  corresponding  to 
significant  structure  rearrangement  and  the  sones  where  the  inner  stress  doesn’t 
depend  up<m  the  duster  extension.  The  spasmodic  transitions  from  one  sone  to 
another  (fig.2)  reflect  the  spasmodic  transfmtmations  of  the  system  of  the  bonds. 

For  aluminum  system  we  investigate  the  influence  of  media  pff  on  its  deformation 
properties.  The  interaction  of  AlM-clnster  with  water  mdecules,  as  well  as  with  the 
additional  prot<ms  or  hydroxyl  groups  is  examined. 

The  dectron-donor  adsmbates  (  like  water  or  hydroxyl)  can  coordinate  m  the  ”on 
•top*  position  on  aluminium  atinns  with  sp’-hybridisidion  of  orbitals  (  like  (111) 
and  (100)  ffices  )  or  in  the  "bridge”  or  the  so-calfed  "slanting  bridge”  position  on 
the  atoms  with  the  proper  sp*-  hybridisation  of  mbitals  (  like  the  (110)  &ce).  In 
order  to  reafise  the  adsorption  of  dectron-acceptor  adsorbate  (like  hydroxoninm)  the 
complete  rearrangoaoent  of  dectroo  density  is  needed. 
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Pennatioii  bMt*  wad  intaratmnie  fore«« 
vcrtnw  bMiie  eliutUr  Al2e  exUiuilen 
with  one-parUcal  adaorbaUa 
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The  comparison  of  the  properties  of  the  systems  containing  1  •  4  water  mdecules 
shows  that  though  the  stability  of  the  clusters  under  investigation  increases  in  the 
presence  of  the  additional  water  molecules;  system  behavior  under  deformation  does 
not  depend  upon  the  number  of  water  molecules  if  it  is  greater  than  3.  The  thing  is 
that  the  optimised  cluster  Al^t  —  nH^O  can  make  chemical  bonds  only  with  3  water 
molecules  while  the  forth  one  is  coordinated  by  the  hydrogen  bond  in  the  second 
sphere. 

The  adding  of  a  proton  to  a  basic  cluster  makes  the  latter  less  stable  (fig.  2),  and  as 
it  occupies  the  ^bridge” -jfonHoa  on  two  aluminium  atoms  the  bond  between  them  is 
broken  and  the  distance  Ali  —  AIu  in  the  free  protonised  cluster  is  a  little  bit  longer 
than  in  the  banc  one  (fig.2).  The  first  aone  of  dastidty  as  well  as  the  decohesion 
energy  is  also  perceptibly  luger. 

In  the  cluster  moddliag  the  acidic  medium  the  HtO*  -ion  does  not  make  a  chemical 
bond  with  any  aluminium  atom  and  the  interaction  is  reduced  to  a  hydrogen-like 
coordination.  The  character  of  the  force-curve  (  fig.3)  is  very  close  to  that  for  the 
basic  cluster  and  in  the  stretching  sone  of  defo>rmati<m  almost  coincides  with  the 
curve  ka  the  duster,  connected  with  2  water  molecules:  so  the  add  medium  seems 
to  make  no  ptindple  influence  on  the  stretching  of  the  duster.  In  the  compression 
sone,  on  the  contrary,  the  work  of  deformation  is  visibly  smaller  than  both  for  bade 
and  for  water-containing  dusters. 
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Interatomic  forces  versus  basic  cluster  Algg  extension 
in  neutral,  acidic  and  alkaline  media 


Pig.  3 

The  alkaline  medium  greatly  influences  the  deformation  properties  of  the  aluminium 
cluster.  Though  the  arid  hydroxyle  group  greatly  increases  the  stability  of  the 
cluster,  the  deformation  properties  are  similar  to  that  of  the  protonixed  cluster 
(fig.2).  The  situation  completely  changes  when  we  model  the  sUghtly-alkaline  medium 
by  adsorbing  on  the  base  aluminium  cluster  a  hydroxyle  ion  together  with  a  water 
molecule.  In  this  case  the  first  sone  of  elastic  deformations  is  spread  on  the  large 
intenral  of  external  tensions  (fig.  3)  and  the  work  of  deformation  is  much  higher  than 
for  all  the  previous  systems. 

So  while  the  acidic  medium  reduces  the  work  necessary  for  the  aluminium  cluster 
destruction,  the  alkaline  medium  expands  the  sone  of  elastic  deformations. 

CONCLUSION 

Here,  the  object  was  to  demonstrate  the  potentialities  of  semiempirical  cluster 
approach  to  quantum-chemical  calculations  of  RE.  The  results  presented  above  should 
be  treated  as  the  examples.  The  semiempirical  cluster  methods  have  an  advantage 
over  the  "Ah  initio"  ones,  because  the  former  allow  to  deal  with  significantly  larger 
systems  containing  up  to  350  atoms.  This  enables  to  calculate  directly  the  structure 
rearrangement  eflbcts  in  the  crack  tip  under  the  influence  of  different  environments 
and  external  stress. 

However,  care  must  be  exercised  in  the  interpretation  of  the  results  obtained  by 
sendeminrical  methods.  It  should  be  recognised  that  the  last  are  no  more  than  rough 
approximation.  Neverthdess,  due  attention  should  be  given  in  our  opinion  to  the 
api^cation  of  semiempirical  cluster  methods  to  RE  investigation. 
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The  repmt  consists  of  three  parts.  First  of  all,  die  dynamics  of  a  single 
macromolecule  in  a  system  of  entangled  weakly-coupled  polymer  chains  as  a 
foundation  of  theory  of  viscoelasticity  for  polymer  solutions  and  melts  are 
considered.  In  die  second  part  consequences  of  maoomolecule  dynamics  for  the 
dieoiy  of  constitutive  equation  for  a  system  in  whole  are  discussed.  In  the  third 
part  method  of  calculation  of  relaxation  properties  of  sur&ce  layer  in  polymer 
composites  is  demonstrated. 

DYNAMICS  OF  MACROMOLECULES 

The  thermal  motion  of  macromolecules  in  the  non-dilute  solutions  and  melts  and 
its  coimection  widi  irreversible  processes  in  systems  of  entangled  polymer  chains 
were  investigated  over  the  last  decades.  Different  research  groups  have  made 
valiant  attempts  to  find  a  unite  approach  to  polymeric  systems  that  would  provide 
an  explanation  and  description  of  nonequilibrium  phenomena  in  entangled 
polymers  fiom  a  unified  point  of  view. 

Essential  progress  in  undmstandiug  the  diennal  motion  of  die  macromolecules  was 
achieved  by  de  Gennes  [1]  and  by  Edwards  {2]  during  die  last  two  decades,  in  the 
last  years  it  was  understo^  diat  the  localization  of  a  macromolecule  in  "a  tube", 
postulated  by  Edwards  and  de  Geimes,  follows  from  a  formal  iqiproach  based  on 
the  most  general  assumptions  in  the  dynamics  of  a  macromolecule  due  to  Russian 
researches.  Apart  of  die  localization  effect  die  approach  spears  to  be  rather  broad 
to  give  explanations  to  all  known  facts  of  dynamics  of  weakly-coupled 
macromolecules  (viscoelasticity,  optical  birefrengence  and  so  on)  that  was 
demonstrated  earlier  (see  review  [3]). 

It  is  convenient  as  die  starting  point  of  the  analysis  to  take  an  approach  which 
appeared  in  the  papers  by  Russian  researches  [3].  According  to  die  Russian  theory 
the  nei^bours  of  the  chosen  macromolecule  form  on  average  a  deformable 
relaxation  medium.  The  chosen  starting  point  is  based  on  the  most  general 
assimqrtions  and  determines  the  first  -  order  approximation  to  the  dynamics  of 
macromolecule. 

We  consider  concentrated  solutions  and  melts  of  linear  polymers  as  systons  of 
weafcfy  connected  macimnolecules:  every  chain  can  change  its  neighbours  and 
moves  among  the  odier  chains.  When  we  consido'  relatively  slow  motion  of 
^stem,  eadi  macrmnolecule  is  schematically  described  as  consisting  of  N+1 


lineally  coupled  Brownian  particles,  so  that  we  shall  be  able  to  look  iqxm  the 
system  as  on  a  suqwnsion  of  many  interactiiig  Brownian  particles  suspend  in  a 
viscous  "mmiomeric  liquid". 

We  consider  the  motion  of  the  described  system  under  deformation  widi  a 
constant  velocity  gradieat  tensor 

vy=(avjy(ar.)  (i) 

so  duit  a  particle  located  at  r“is  dragged  with  mean  velocity 
A  force  acting  on  any  part  jf  the  system  gives  rise  to  the  excrmon  of  ^e  entire 
ensemble  of  Brownian  particles,  so  that  when  behavior  of  the  system  is 
investigated  we  have  to  consider  die  collective  motion  of  all  the  particles  in  the 
same  way  that,  for  example,  we  examine  the  motion  of  the  ensemble  of  atoms  in 
solids.  Our  task  is  therefore  to  find  the  normal  coordinates  of  the  polymer  system, 
i.e.,  die  variables  that  vary  independently  of  one  anodier. 

The  identification  of  the  normal  coordinates  can  be  carried  in  two  stages,  bearing 
in  mind  the  particular  properties  of  the  system  (strong  interaction  along  the  chain 
and  weak  interaction  between  die  macromolecules).  The  task  of  the  first  stage  is  to 
determine  the  dynamics  of  a  single  macromolecule,  surrounded  by  all  the  others. 
To  implement  the  first  stage,  we  must  eliminate  all  variables,  other  than  diose  diat 
refer  to  the  cnosen  macromolecule.  This  procedure  is  not  too  simple,  but  we  can 
imagine  the  general  form  of  the  fins.'  results  before  the  calculations  are  carried  out. 
The  requirements  of  proper  covariance  and  of  linearity  in  coordinates  and 
velocities  determine  [3]  the  general  form  of  the  equation  for  the  dynamics  of  the 
chosen  macromolecule. 

This  is  followed  by  the  second  stage  m  which  the  normal  coordinates  of  an 
individual  macromolecule  are  determined  [3]. 

The  implementation  of  the  program  gives  stochastic  equations  for  difiusion  and 
relaxation  modes,  respectively 


m{<Pp^{dP)  =  -  Jp(s)(p  “  -  Vj,p ,“)  ds  +  o^{t)  (2) 

m{(P pPy{d?)  =  -  fP(s)(Pj“  -  VyPj“)j_g  ds  - 

-  WXPitt  ^  -  ^P^aP”  +  '^“(0. 

a=1.2....,N  (3) 

where  m  is  die  mass  of  a  Brownian  particle  associated  with  a  piece  of  the 
macromolecule  of  length  M/N,  2\xT  is  elasticity  coefficient,  T  is  tenqierature,  oa 
Q=//2(vy-V|.).  For  large  Nani  small  a  the  eigenvalues  are  then  given  by 

X  =  (naJN)^, 


a=l,2,...,«N 


(4) 


The  fint  teams  on  die  right  of  die  Eqs.(2)  and  (3)  ne  dw  hydroefynamic  dn^ 
fioeees,  Ae  second  tenn  m  Eq.(3)  fqmsents  die  mtmndeciilar  leaslaace  (be  to 
the  diange  in  die  sluqie  of  die  macnmioleinilar  coil  (kinetic  stiQhess)  the  diird 
teim  in  Eq.(3)  represent  die  elastic  fence  dne  to  nearest-neighbour  Brownian 
particles  along  the  chain.  The  last  terms  in  Eqs.(3)  and  (4)  are  die  random  tfaennal 
forces  whose  statistical  pnqierties  are,  as  amal,  defin^  so  that  the  eepulilnium 
values  of  die  calculated  quantities  are  die  same  as  diose  already  known. 

According  to  [4]  we  should  also  write  down  dw  drag  and  resistance  forces  in 
Eqs.(2)  and  (3)  in  dw  form  of  nonlocal  eiqnssions,  since  exitaticMi  directly 
duou^  dw  chW  propagates  to  a  distance  <k>,  i.e.  a  distance  that  is  large  in 
conqiarison  widi  die  size  of  dw  Brownian  particle  under  consideration.  However, 
for  the  sake  of  sinqilicity,  diis  will  not  be  done  here,  aldiough  die  consecpiences  of 
a  nonlocal  effect  be  noted  latter.  Eipiations  (2)  and  (3)  also  rely  on  the 
questionable  assumpdim  dud  each  particle  in  the  chain  is  in  die  same  isotropic 
situation. 

Apart  of  these  notes,  Eqs.(2)  and  (3)  are  the  most  general  equations  for  the  low- 
fiequency  modes  of  dw  system  in  dw  case  of  a  linear  dependmee  on  variables. 
Howevtr,  dw  memory  functions  P(s)  and  (p(j)  can  not  be  determined  from  general 
considerations:  dwy  must  be  found  by  calculations  such  as,  for  example,  those 
reported  in  [5],  or  sinqile  heuristic  consideration  must  be  abandoned  as  was 
observed  in  [3]. 

In  a  simple  case,  die  one-sided  Fourier  transforms  of  the  memory  functions  were 
found 


=  +B/(l-/<at)),  (p[sl  =  C£/(7-«ox)  (5) 

where  C  is  the  monomer  friction,  B  is  the  measure  of  increasing  of  friction  due  to 
neighbouring  macromolecules,  E  is  the  measure  of  "internal  viscosity”  due  to 
entanglements  of  macromolecules  in  the  system.  The  relaxation  time  x  in  (5)  is  the 
relaxation  time  of  the  matrix  which  we  later  calculate  as  the  main  viscoelasticity 
relaxation  time  of  the  entire  system  in  a  self-consistent  manner. 

The  parameters  >^JE  and  x  are  the  phenomenological  parameters  in  the  theory. 

We  can  find  B~K^ ,  where  M  is  die  lengdi  of  macromolecule,  if  we  take  the  non¬ 
locality  of  medium  into  account.  According  to  [4]  5=2,  experiments  give  a  value 
of  die  exponent  slighdy  more  than  2,  consistent  with  model. 

Equations  (2)  and  (3)  can  be  looked  iqxm  as  the  first-order  qiproximation  to 
nonequilibrium  phenennena  in  undihite  polymers.  This  excludes  effects  due  to 
nonlinear  terms,  e.g.,  dw  reptadon  phenomena  noted  by  de  Gmmes  [1],  which 
qqwar  when  die  difference  between  mobilities  along  and  at  right  an^es  to  die 
(diain  is  taken  into  account  This  can  be  described  by  terms  of  order  greater  than 
dw  first  in  a  rigorous  dieory  of  motion  of  a  macrrmudecule  in  a  concentrated 
tystem.  However,  before  we  tom  to  this  problem,  we  must  mumune  the  effects 
represoited  Ity  lirwar  terms. 
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MOLECULAR  NON-LINEAR  THEORY  OF  VISCOELASTICITY 


The  described  results  form  a  basis  on  which  tbe  rheological  behavior  of  polymer 
liquids  can  be  considered.  Now  we  are  in  position  to  formulate  crmstitutive 
equations  for  po^mer  liquids.  In  &ct  we  have  a  set  of  constitutive  equations  due 
to  various  rqiproximations. 

The  equatirms  of  motirm  for  a  set  of  Browning  particles  suqwnded  in  a 
viscoelastic  liquid  are  fmmulated,  which  is  equivalent  to  die  derivation  of  the 
expression  for  the  stress  tensor  which  we  write  in  terms  of  modes  (3)  in  die 
foUowing  form  [3] 

ajk(^>=  <P  “(0Pk“(0>  -  + 

+  (7/7)  J<p(s)[<Pj“(r)pj^Vs)>  <Pj“(t-s)Pj“(/)>ldr}  (6) 

where  p  is  the  pressure,  n  is  the  macromolecular  number  density. 

The  angle  brackets  represent  averaging  over  the  ensemble  of  realization  of  random 
forces  in  the  mode  equations  (3). 

It  is  convenient  to  introduce  internal  relaxate  variables  and  “  to  rewhrite 
stress  tensor  in  another  form 


The  first  set  of  variables  is  connected  with  deformation  of  macromolecule  in  the 
system.  The  second  set  of  variables  reflects  stresses  under  deformation  of  "die 
tube".  It  was  shown  diat  two  sets  of  relaxate  variables  give  the  adequate  picture  of 
relaxation  of  polymer  systems  [3], 


It  was  shown  [6,7]  that  relaxation  equations  for  x.  -^  and  follow  firom  mode 
equation  (3) 


("'“V  -  (*«■' "  'kj  ■ 


-  -  ‘\l%"  ■  Hi  *  (-kj“  -  Hi 


(OOO./  -  -  -  (Ih  1 

*  ■  HS  8^/  -  Wik“> 


41 


Aoxiliaiy  quaitthies  b^,  c^,  “  and/-.  “  are  defined  [7]  as  fiinctions  of 


parameters  of  die  diedfy  ana  anisovqyjriensor 


Relaxati<m  equatums  (8)  and  (9)  cmitain  relaxation  times 

T,  x^  =  x/2  +  x^(B+E),  x^  =  x*/a^. 

It  was  shown  [3,6,7]  that  the  syston  of  constitutive  equations  (7),  (8)  and  (9) 
describes  various  no-linear  effects,  in  particular  anisotropy  of  relaxation. 

The  constitutive  equation  takes  a  simpler  form  for  fluid  polymer  system  [7], 

MACROMOLECULE  IN  A  SURFACE  LAYERS 

The  properties  of  polymer  near  die  surface  of  solid  is  quite  different  fi'om 
properties  of  bulk  polymer.  It  is  necessary  to  take  into  account  this  effect  under 
the  constructing  of  polymer  solid  composite  widi  large  content  of  solid  phase.  In 
the  last  case  almost  all  amount  of  poljmer  is  in  modified  state  in  surface  layers. 
We  can  use  the  described  theory  of  viscoelasticity  to  determine  the  alternation  of 
polymer  properties  in  a  surface  layer. 

As  a  starting  point  we  take  the  dynamics  of  macromolecule  (2)  and  (3),  whereby 
the  interaction  of  particles  with  surface  is  taken  into  account.  The  modified 
equation  of  dynamics  is  non-linear  now  and  we  have  to  use  computer  simulation 
methods  to  came  to  results.  Computer  mediods  allows  us  to  follow  the  stochastic 
trajectory  of  every  Brownian  particle  (bead)  of  chain. 

The  data  obtained  fiom  corrqniter  simulation  are  used  for  calculation  of  the 
correlation  functions  and  mean  square  displacement  of  separate  beads  in  die 
macromolecule.  The  results  should  yield  difteion  coefficient  and  relaxation  times 
of  die  macromolecule. 
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INTRODUCTION 

Conq)uter  simulation  is  widely  used  for  studying  die  mechanical  properties  of 
polymers  because  of  dieir  practical  applications  and  there  is  a  strong  need  for 
accurate  prediction  of  properties  without  resort  to  expensive  experiments. 

It  is  well-known  three  meAods  in  die  modeling  of  polymer  materials.  They  are  die 
stochastic  dynamics  methods  (SD).  The  Brownian  dynamics  mediods  (BD)  and  the 
molecular  dynamics  methods  (MD). 

The  most  accurate  treatment  in  principle,  ought  to  he  MD  method,  because  of  the 
ex{dicit  description  of  atomic  motions;  it  takes  into  account  the  interactions 
between  polymers  chains,  die  interaction  with  solvent  particles  and  surfaces.  It  is 
very  impotent  for  investigation  of  polymer  composite  materials. 

Our  aim  is  to  fnesent  the  original  molecular  dynamics  program  MOLSIM  and 
discuss  their  applications. 

DETAILS  OF  ALGORITHM 

The  technique  of  molecular  dynamics  simulations  is  well-known.  The  classical 
equatitms  of  moticm  are  solved  for  a  system  of  interacted  particles.  Basic 
[Hincqiles  of  the  method  have  been  described  by  Raman,  Parrinello  (constant- 
pressure  technique).  Nose,  Hoover,  Andersen  (constant  temperature  algorithm). 
Then  tins  method  have  been  develcqied  by  Berendsen  (loose-coiqiling,  constant- 
inessure  technique)  and  extended  by  Brown  and  Cluke  for  use  on  polymer 
samples. 

This  ideas  form  the  basis  of  the  MOLSIM  program. 

The  most  inqiortant  features  if  die  MOLSIM  program  are  original  algorithm  of 
searching  neighbors  for  evaluating  nonbonded  interactions  and  database  FIELDS 
of  empirical  parameters  of  die  Force  Field  of  interatomic  interactions. 

Original  algorithm  of  effective  regional  search  of  neighbors  is  realized  in 
MOLDYN  module  of  MOLSIM  program.  Operational  feature  of  program  are  the 
core  memory  requirement  depends  only  on  die  number  of  particles  and  not 
depends  on  die  value  of  cutoff  radius.  This  makes  it  possible  to  simulate  efficiendy 
very  large  system  and  use  long-range  potential. 

The  database  pn^jam  FIELDS  provide  to  formation  a  force  field  describing 
interactions  between  atoms  (particles)  generally  contains  term  of  covalent  bond 
interactions,  btmd-angle  interactions,  torsion  angle  interactions,  and  nonbonded 
intmactions.  The  database  contains  variable  potential  fimctions.  And  widely  range 
of  potential  fimction  are  also  realized  in  MOLDYN  inogram. 


44 


MCXJSIM  PRCXatAM  EffiSOUPTION 


Prognun  NKX.SIM  haa  beta  developed  to  provide  bilenclioii  Foroe  Fields 
Fonnirtion  for  given  moleciilsr  structures  and  molecular  tlynaimcs  shnalrtions. 
This  program  ctmsisis  of  SUBGR,  FIELDS  and  MOLDYN  modules.  Functional 
sdmne  of  program: 

File  •.MOL  -*  SUBGR  ->  FIELDS  -►  MOLDYN. 

SUBGR  module  is  used  fot  analyzing  given  molecular  stnictures  and 
determination  all  isomondiic  subgraidis  (uncolored)  of  certain  tautologies  and  dw 
colmed  subgrqths  of  the  same  tc^logies  in  a  colored  grqth,  which  represent  one-, 
two-,  tree-  a^  four-botfy  atannic  chains.  "Color”  is  connected  wiA  an  atomic 
name  or  an  atannic  groqt  name. 

FIELDS  module  provides  access  to  the  Data  Base  of  empirical  interatomic 
potential  fields  to  allocate  the  fmce  fields  parameters  for  the  constants  in  the 
potential  energy  fimction  for  given  molecular  structure.  It  nuQ'  be  used  also  to 
create  a  new  potential  fields  for  users.  Database  provides  to  create  a  Force  Field 
vdiich  involves  valence  and  nonvalence  intnactions 


where  E  .  includes  the  terms  arising  fiom  covalent  bond  fcmnation  and  E^^^ 
includes  tte  longrange  noncovalent  interactions. 

As  usual,  E^  includes  Umd  stretdi,  angle  bend,  dihedral  angle  torsion,  one-center 
angle-angle  cross  terms,  v^e  die  nonbond  terms  consist  of  van  der  Waals  and 
electrostatic  terms. 

The  MOLDYN  {Hogram  has  been  developed  for  molecular  dynamics  simulations. 
The  conqniter  simulation  techniques  based  on  the  well-known  molecular  dynamics 
(  MD  )  method  and  can  be  used  fmr  studying  the  thermocfynamical,  structural  and 
(fynamical  properties  of  polymm'  matoids.  The  program  MOLDYN  solves  the 
classical  equatimis  of  modon  of  a  3D-system  of  mutually  interacting  particles. 

The  equations  of  modrni  is  solved  by  using  of  the  Veriet  algorithm.  The  desired 
tenqrenture  is  controlled  via  scaling  periodically  die  velodties  of  the  atoms. 

The  program  MOLDYN  calculates  die  thermodynamic  quantities:  kinetic, 
potent  and  total  system  energy,  temperature  ant  oAer  medumical  values  of  die 
systm. 

Progtanuning  is  FORTRAN  77.  The  computations  may  be  pofmmed  on 

die  VAX  hardware  as  well  as  on  IBM  PC  AT/XT  computers. 

SIMULATIW 

KRdecular  dynamics  conqmter  simulatioos  have  been  used  to  stu^  die  extantional 
(Young's)  modulus  of  oriented  linear  ptdymer  chains  nd  poisoner  samples 
generated  by  uang  a  modified  sdf-avoiding  random  walk  techniqne.  Investigations 
were  pctfimned  under  constant  tenskm  at  a  range  of  tcnqierature. 
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The  detuled  tiwh—iffwi  of  the  load  dynamics  of  pcdymer  chains  and  die  surface 
bound  polynier  chains  have  been  investigated. 

hi  die  sidbse(|Beat  we  aasuaie  to  extend  ow  simulatum  <»  polymer  con^)osite 
materiab  investigations. 
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Thoe  we  two  types  of  the  periodical  systems,  one  of  Aem  is  die  based  on  the 
properties  of  electrons,  diat  were  determined  by  periodic  system,  and  another  •  on 
dm  nucleus  properties.  Periodic  system  is  one  jvojecdon  of  die  whole  periodic 
properties  a^  another  projecdrm  is  periodical  scale  of  die  nucleus  properties. 
Many  fimts,  which  can  not  be  described  by  periodic  system,  may  be  seen  in  the 
natural  systems  (i.e.  in  geochemical  structures). 

On  this  base  we  realise  new  directions  in  the  calculations  of  composites  and 
quasinatnral  stractwes  for  the  qiplications  in  electronics,  hi^  temperature 
siqpefcooductor,  materials  w^  extre^  qualities. 
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The  design  and  febricadon  of  |»ef(»ms  for  advanced  composites  have  gained 
ctmsiderable  attention  in  the  lig|it  of  die  recendy  develop^  diree-dimensional 
textile  inefonning  techniques.  Three-diinensional  textile  structural  composites, 
based  upon  braiding,  weaving,  knitting  and  stitching,  provide  enhanced  stiffiiess 
and  strength  in  the  diickness  diiection,ini{»roved  damage  tolerance,  and  most 
importandy,  the  potential  of  near-net-dii^  forming  of  components  widi  complex 
surface  contours.  The  ability  in  die  integrated  and  system  qiproach  to  textile 
structural  composites  from  preform  processing/manufacturing,  matrix 
consolidation  to  composite  charactoization,  and  analytical  modeling  of  composite 
performance  has  greatly  broadened  die  potential  of  textile  composites  in  primary 
structural  triplications. 

Two  aspects  of  the  computer  syndiesis  effort  in  die  modeling  of  fiber  architecture 
and  performance  of  textile  composites  are  outlined  below.  First  in  the  area  of 
process  simulation,  die  develoinnent  of  a  computer  simulation  of  a  general 
Cartesian  braiding  process  has  led  to  an  understanding  of  process  fimdamentals. 
The  simulation  allows  for  the  identification  of  individual  yam  paths,  number  and 
location  of  yam  groiqis,  and  braid  geomedy.  Innovative  braid  geometries  were 
simulated  aib  the  prefmms  fabricated  to  demonstrate  the  feasibility  of  fabricating 
a  wide  range  of  fuefonn  architectures  given  an  advanced  braiding  machine. 
Addidoiudly,  interestii^  distributions  of  yam  groiqis  have  been  found  which 
Sliest  an  iqqilicatitm  to  l^lnid  ctmqiosites.  These  new  and  advanced  braids, 
termed  "multi-step"  braids,  are  only  possible  wife  individual  row/column  control. 
The  multi-step  braiding  process  gr»fnty  extends  fee  range  of  possible  preform 
microstractures. 

Next  fee  feennal  and  mechanical  property  modeling  and  prediction  follows  a 
hierarchical  qqnoach  which  entails  meclu^cal  and  geometric  modeling  of  each 
fiber  tow,  followed  by  fee  identification  of  mechanical  relationships  between  each 
of  fee  tows,  and  fee  use  of  an  assembly  scheme  to  determine  the  feermo-elastic 
properties  of  fee  contyosite.  The  predict  stiffoesses  for  typical  30  architectures 
are  dun  ctmqiared  wife  fee  experimental  data  to  validate  fee  model  predictions. 
Estimates  fin  fee  stroigfe  of  fee  textile  ccnityosites  along  fee  princi^  material 
directions  are  obtained  from  fee  tow  stiengfe,  fee  fiber  volume  distribution,  and 
fee  parameters  characterizing  fee  tow  undulation. 
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Throu^  the  analytical  simulaticns  and  detailed  parametric  studies,  a  systematic 
iq)proach  has  be»  accmnpindied  to  understand  die  relaritmships  between  key 
mataial  and  gemnetric  parameters  and  die  resultant  composite  properties. 
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INTRODUCTION 

The  study  of  nonstationaty  wave  process  in  explosion  or  shock  is  of  considerable 
intoest  for  a  multiplicity  of  uses.  As  a  rule,  the  solving  of  such  problems  is  associated 
with  considerable  mathematical  and  computing  difficulties.  This  is  particularly  true 
in  regard  to  the  media  with  complex  behavior  under  the  high-speed  loading.  However, 
there  exists  class  of  problems  when  difficulties  arise  even  in  the  case  of  simple  media 
such  as  ideal  gase. 

Asymptotic  methods  are  very  useful  in  qualitative  and  qusmtitative  analysis  of 
complex  wave  processes.These  methods  are  based  on  expsmsion  of  solution  in 
asymptotic  series  over  pven  smaU  parameters  set.  However,  direct  expansion 
never  leads  to  desired  result,  because  of  wave  problems  specificity.  The  reason 
is  that,  the  validity  domain  of  direct  expansion  in  space  of  independent  variables 
sqtpears  to  be  too  small.  The  overcoming  of  this  difficulty  is  possible  by  using  of 
many-scales-method  (MSM).This  method  was  developed  in  1970-s  by  the  efforts  of 
many  investigators  (Cote,  Kevorkian,  Lick,  Taniuti  and  Wei,  and  othere).  The 
contemporary  form  of  MSM  supposedly  was  suggested  by  Leibovich  and  Seebass, 
and  Oikawa  and  Yajima. 

In  nonlinear  wave  theory,  MSM  is  used  to  decompose  initial  complex  system  of 
equations.  As  a  result,  the  independent  system  of  equations  for  the  functions,  which 
are  Riemann’s  invariants  analogs  has  been  derived.  In  first-order  (nonuniformly 
vaHd)  iq;>proximation  these  functions  are  constant  along  characteristic  directions. 
Keefnng  an  account  of  nonfinearity  and  kinetic  processes  leads  to  their  slow  variation. 
The  necessary  condition  of  MSM  validity  is  weak  nonlinearity,  that  is,  inrignificant 
variaiton  oi  substance  denaty  in  a  wave  or  small  ratio  of  pressure  to  bulk  modulus. 
For  the  representative  problems  on  impact-  or  explonon-  induced  deformation  of 
metals,  ceramics  or  composites,  that  condition  is  not  too  restrictive. 

In  the  problems  on  strong  blast  wave  propogation  in  gases  or  in  porous  media 
(foe  example,  in  rocks),  the  rituation  is  fundamentally  new.  At  the  initial  stage 
of  explosion  in  gases,  the  pressure  in  a  shock  wave  extremely  exceeds  the  initial 
pressure  in  a  medium.  In  a  porous  media,  the  pressure  is  extremely  higher  than 
the  compresnon  strength  limit,  which  leads  to  irreversible  paddng.  In  both  cases, 
the  nonfinearity  is  strong  and  MSM  can  not  be  used  in  the  general  formulation. 
The  most  effective  method  for  qualitative  study  of  this  type  of  problenu  is  the 
representation  of  these  problems  as  of  models  with  lumped  parameters  by  using 
of  asymptotic  expanrions.  The  fundamental  results  in  this  field  have  been  obtained 
in  1960-60s.  We  can  state  with  assurance  that  modds  with  lumped  parameters 
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lum  kwt  popnhiitT  <in«  to  extenacm  d  ocunpiiting  potcntiaHtiet.  However,  for  the 
optiaul  eoet^  problane,  theae  modele  appear  to  be  the  only  avilable  method  of  the 
■idMeaHj  Ml  qoMtative  etadj. 

Here,  we  the  poeailMfitiee  aad  pecnHaiitiet  of  application  of  MSM  and 

modeli  with  himped  parametera  to  itndping  nmutatfonarj  wave  ptoceaaea  under 
ahock  and  explockm.  The  «TaTwpW  of  qualitative  inveatigati<m  and  d  solving  of 
tidied  proUons  are  discasaed. 

GENERAL  PROCEDURE  OF  DECOMPOSITION 

A  host  of  models  of  continuous  media  results  in  system  of  equations  of  the  type: 

^  i  =  1 . N  (1) 

were  u  is  the  unknown  vector-function,  Cik^S)  is  NxN  matrix,  D4u]  -  are  certain 
operators,  p  «  1  is  a  small  parameter.  In  particular,  system  (1)  describes  one¬ 
dimensional  nonlinear  waves  in  visco-elastic  and  elastoplastic  media  with  intrinsic 
degree  of  freedom,  in  relaxating  fluids  and  gases,  membranes,  etc. 

The  proUems  on  shock  and  explosion-induced  deformation  of  materials  are  usually 
stated  as  the  boundary-value  problems.  However,  it  is  more  convinient  to  consider 
the  general  procedure  of  decompoation  on  the  example  of  Canchy  problem.  This 
allows  to  demonstrate  that  procedure  in  the  most  general  form,  while  the  extension 
of  the  result  to  the  boundary-value  problems  presents  no  difficulties. 

Let  if  s  «(0,  x),  -00  <  X  <  00,  be  given  at  the  initial  moment  t=0,  with 

||»(0,»)-ttc||'ve<l,  (2) 

where  ««  is  unperturbed  homogeneous  stationary  state  of  the  medium.  It  is 
reasonable  to  represent  the  solution  of  system  (1)  in  the  form  of  expansion  near 
the  homogeneous  stationary  state 


u  =  ttc  +  eui  +  «*uj,  +  + ... 


(3) 


.  It  is  supposed  about  Cik(it)  and  that 


+ ... 


Dilu)  =  eD^‘^(ue,wi]  + ... 


(4) 


By  substituting  the  expansion  (3),(4)  in  (1)  we  find  the  system  of  equations  in  uk. 

The  first  stq>  of  decomposition  procedure  involves  diagonalisation  d  that  system. 
The  possilnlity  of  diagonalisation  is  provided  by  the  reality  of  the  eigenvalues  of 
matrix  C^\«c).  Tb  put  it  otherwise,  to  a  first  iq>praximation,  system  (1)  is  supposed 
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to  b«  a  hyp<**»*^t*  ^  ***  mtioduee  the  nniUry  matrix  Tji{Se)  and  vector  Vj 

■o  that  TiiC^^Tu  =  Tum  when  Xji'Sc)  an  eigenvalues  o£  matrix 

Ca(d^).  expansion  of  Vj  over  small  parameters  foQows  £tom  (3),  wen  all  Vu 
an  alosig  thsir  eharactsristk  directioas.  At  the  second  step  it  should  be 

adopted  that  aB  Vi^  along  their  diaracteristks  an  slow  functions  of  time,  that  is 
Vif  has  a  form  Vy  =«  with  =  x  -  Xjt,T,  =  et,T^  =  (it.  Then,  the 

equations  for  the  seoondswto  terms  of  expansion  an  the  following: 


avij 

5t. 


svif^ 

dl 


=  TjiD\ 


h(^hv. 


(Vil 


dVij 

dr^ 


(5) 


The  derivatives  in  the  left  parts  of  equations  (5)  are  calculated  with  fixed  (j,T„T^ 
while  in  the  right  parts  Vn  =  Vifc((j  +  AAjst,r«,r„),  AA^s  =  Aj  -  A*  . 

FVirthermom,  when  integrating  (5),  we  must  eliminate  the  terms  which  increase  with 
the  "fast”  time  t.  This  can  be  arranged  by  the  selection  of  "slow”  time-dependence  (i. 
e.,  r«,r^  —  dependence  )  of  functions  Vij.  Thenby,  the  smallness  of  the  second-order 
terms  in  comparison  with  first-order  ones  in  expansion  (3)  can  be  extended  to  the 
asymptotically  long  time  interval  i  <  0(e~‘  -t- 

For  the  application  of  interest  here,  the  most  frequently  used  type  of  system  (1)  is 
that  of  real  one,  where  are  linear  comlnnations  of  uj  components  and 

are  linear  operators.  We  assume  these  operators  to  preserve  the  order  or 
function  value. 

Feasibility  of  decomposition  is  determined  by  the  character  of  phase  velocities 
spectrum  {Aj(t(e)}.  Let  us  suppose  that  there  are  no  closdy-spaced  eigenvalues  in 
the  spectrum  of  phase  velocities;  that  is  the  condition  |AA;k|  >>  max{ey(i,}  should 
be  met  for  any  j  ^  k.  Uniformity  requirement  for  the  first  terms  of  expansion  (3) 
can  be  ssdisfied  if  the  initial  disturbance  is  restricted  ,  for  example,  by  the  condition 

I  /  -  «*i]|  ^  *  (®) 

for  any  xi ,  Xj  €  (—00,00).  As  the  consequence  of  (5)  and  (6),  the  following  estimate  r 
can  be  obtained: 

Here  Qi  are  nonnegative  constants,  Qi  «  In  this  case  the  general 

form  of  the  equations  for  the  first-order  (^prcodmation  (uniformly  valid  at  the 
asymptotically  large  time  interval)  to  the  solution  is 

^  j  =  1, ...,  W  (7) 
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Here  the  derivatives  with  respect  to  time  t  are  calculated  vrith  fixed  phase  variables 
The  term  is  the  right  side  of  equatioo  (7)  is  the  sununaad  of  linear  combination 
which  does  not  depend  on  the  "fitsf  time. 

The  extension  of  outlined  procedure  to  the  case  of  closely-spaced  eigenvalues  presence 
is  not  difilcult.  Let  only  two  eigenvalues  As  and  are  close  to  each  other.  Let 
us  introduce  one  further  "slow”  time  rx  =  AAsmt-  Then  +  r*,T„T^)  and 

secular  terms  elimination  leads  to  the  coupling  equations  for  Vn,  and  Vim-  The 
other  (N-2)  equations  of  the  form  (7)  remain  unchanged. 

Modification  to  the  case  of  boundary-value  problem  also  leads  to  the  equations 
of  the  form  (7).  The  minor  distinction  consists  in  using  the  transformation  (  — » 
X  and  (j  —*  —(j  in  (7). 

The  consideration  carried  out  here  refers  only  to  unidimensional  proUems.  However, 
there  is  a  host  of  interesting  multi-dimensional  problems.  The  impact-induced 
deformation  wave  propagation  in  solids  is  an  example.  Within  the  limits  of  MSM,  the 
natural  way  to  consider  multi-dimensional  effects  is  the  using  of  so-called  Leontovich 
pwabolic  approximation.  To  our  knowledge  this  approximation  was  first  applied  to 
nonlinear  problems  by  Zabolotskaya  and  Hohlov,  and  Kadomtsev  and  Petviashvili. 

DEFORMATION  WAVES  IN  THE  ELASTOPLASTIC  MEDIA  : 
"INTERNAL”  SMALL  PARAMETERS 

Elastoplastic  behavior  is  a  characteristic  one  for  solids  in  impact  and  explosion 
-  induced  deformation.  In  the  range  of  insignificant  shear  stresses,  the  material 
undergoes  elastic  deformation,  while  undex-  the  exceeding  of  certain  threshold  stresses 
the  material  fiows.  In  such  a  situation,  the  application  of  MSM  to  description  of 
deformation  waves  has  certain  peculiarities.  These  peculiarities  are  caused  by  strong 
dependence  of  medium  kinetic  characteristics  on  stresses  (and  temperature)  in  the 
domain  of  elastic-to-plastic  transition.  This  is  a  handicap  to  the  expansion  of  the 
kinetic  functions  into  the  series  over  the  small  parameter  e. 

Let  us  accept  the  initial  unidimensional  equations  of  Maxwell  elastoplastidty  in  the 
form  proposed  by  Godunov  and  Romensldi.  Deformation  kinetics  in  this  model  is 
determined  by  a  single  scalar  function  of  stress  tensor  invariants,  namely,  by  the 
shear  stresses  relaxation  time. 

At  first  the  initial  system  of  equation  must  be  represented  in  the  dimensionless  form. 
The  resulting  set  of  small  parameters  may  be  written  as; 

_  wj  ^  _  3(  +  Arf 

*  ~  Ko"'  Pf,'  3KoTo 

K  -cj 

Here  K9  and  po  are  respectively  bulk  modulus  and  solid  denrity,  and  k 
are  internal  friction  and  heat  transfer  coefficients,  ci  is  phase  velody  of  lon^tudinal 
elastic  waves,  co  =  y/Kt/pt,  <rj  and  re  are  characteristic  values  of  stress  and 
time  from  the  boundary  condition. 
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One  further  small  parameter  S  defines  the  width  of  dastic-to-plastic  transition  domain 
in  the  rdation  between  rdaxation  time  and  dominant  shear  stress.  The  substitution 
of  new  small  parameter  S  provides  the  agreement  between  the  calculated  deformation 
wave  structure  and  its  experimentally  observed  splitting  into  elastic  and  plastic  waves. 
It  is  i^iparent  that  the  expansion  <d  the  rdaxation  equations  over  the  t/S  ratio  makes 
no  sense  when  f  <  e  .  Hence,  the  initial  form  of  rdaxation  equations  must  be 
preserved.  The  decomposition  procedure  therewith  reproduces  the  discussed  above 
one,  exept  that  the  new  restriction  is  substituded  to  provide  the  validity  at  the  lenth 
range 

It  is  of  interest  to  compare  the  results  of  numerical  solution  of  initial  dastoplastic 
equations  and  M SM  ones.  Such  comparision  was  performed  for  the  explosion-induced 
deformation  wave  in  iron.  The  results  of  numerical  calculations  was  found  to  be  in  a 
good  agreement. 

IMPULSE  PROCESSES  IN  THE  ELASTIC  LIQUID-FILLED  PIPES 


Additional  examples  of  MSM  efficiency  are  the  nonlinear  problems  on  wave  processes 
in  a  simple  hydro-elastic  systems.  Here  we  shall  consider  nonlinear  waves  induced  by 
a  local  build  up  of  pressure  in  the  liquid-filled  elastic  pipe.  Liquid  is  presumed  to  be 
an  ideal  compressible  one. 


On  the  assumption  that  the  shell  deformations  are  insignificant,  the  hydraulic 
approximation  for  the  liquid  in  motion  can  be  used.  The  equation  for  the  relative 
variation  of  pipe  cross  section  under  the  effect  of  angle-independent  loading  may  be 
written  as  .  ,  . 

a>s  E.A»  a*S 


12p«(l  - 1^)  dz*  p,w 


:S+- 


P»ah 


Here  E,fPt,v  •  are  Young  modulus,  density  and  Poisson  ratio  of  shell  material 
respectively,  a  is  radius  and  h  is  thickness  of  the  shell. 


The  possibility  of  MSM  application  is  based  on  the  presence  of  two  small  parameters. 
The  first  is  the  ratio  of  characteristic  pressure  to  the  bulk  modul'os  of  liquid  e  = 
and  the  second  is  liquid-shell  interaction  parameter  p  =  2poCoa/E«h. 

Decomposition  of  the  initial  system  results  in  two  independent  Witham  -  type 
equations  for  the  Riemann’s  invariants  V±  =  («  ±  Ap/poCe)/coc- 


(9) 


where  t  =  4ewot/(l  +  n),(+  =  wt(i  T  iiJco),V9  =  =  p/(l  -(-  n)e,n  -  is 

politropic  exponent  for  liquid  or  gas.  In  equation  (9),  the  integral  operator  kernel  is 
defined  by  its  Fourier  tranCmirm 


&X 


R{w)  = 


1 

l—w^+  D^v* 


(10) 


D*== 


12(l-,^)rfa»c} 


54 


The  eqiuttton  (9),  in  tww  of  (10),  me  ctmaenmUTe  and  can  be  &Nmd  from  the  weO* 
blown  eariatinnal  ptmci^  tat  Wtham’a  equationa.  It  foUowa  that  eqnatkm  (9) 
haa  thtee  potynomial  coaaerratkm  laws  and  the  cmuervatkm  law,  which  cmitaina  the 
independent  variaUea  in  an  explicit  form. 

In  the  hmg'Wave  limit,  eqnatkm  (10)  nuqr  be  leananged  into  the  weQ-atndied 
KdV-cquatkm.  The  nanal  acenaao  at  initial  diatnxbance  endntieo  Car  non-local 
generaliiationa  of  KdV-eqnatiom  it  diatnrbance  decmnpodtion  into  the  sditon  -  like 
Kdntiona  and  damping  background.  In  this  connection,  the  localixed  steady-ctate 
■obtiona  of  equation  (9)  are  of  impmtance.  Theae  aolutiona  depend  only  on  variable 
y  =  (  -1-  oAr,  where  A  ia  aome  arbitrary  conatant. 

In  D  — »  0  limit,  the  atationary  aolutiona  can  be  found  without  difficultiea.  It  tuma 
out  that  localioed  oohton-like  aolutiona  exiat  only  in  the  range  0  <  A  <  1;  there  exist 
two  different  Idnda  of  localised  solationa:  the  omooth  onea  and  the  oingular  onea.  The 
smooth  solutiona  (acditona)  are  characterised  by  poritive  pressure  in  pulse  while  the 
singular  ones  (we  call  them  cavitons)  differ  in  that  they  have  a  negative  pulse  pressure 
and  two  pmnta  with  infinite  derivatives.  The  computer  simulation  has  showed  that 
(with  D  =  0)  the  initial  disturbance  evolution  really  shows  up  as  decomposition  into 
sequence  of  solitons  and  cavitons  and  fast-damping  background. 

With  D  ji  0,  the  localised  solutions  have  been  studied  both  by  numerical  and 
asymptotical  (with  D  «  1)  methods.  The  main  result  is  that  the  spectrum  of 
eigenvalues  A  for  solitons  and  cavitons  becomes  discrete  and  consisits  of  the  finite 
number  of  points  in  the  range  from  0  to  1.  With  small  D,  "the  quantisation  rules” 
are  obtained  fox  the  spectrum  calculation  in  the  first  order  in  D. 

The  continuous  spectrum  decomposition  is  accompanied  by  the  appearance  of 
combined  states  and  weak  stationary  chaos.  The  combined  states  may  be  treated 
as  "molecules”  composed  of  finite  number  of  solitons  and  cavitons.  Stationary  chaos 
is  defined  as  a  series  of  randomly  positioned  sditons  and  cavitons.  Thus  we  call  it 
"soHton  glass”. 

In  adddition,  with  17/0  there  appears  a  new  set  of  localiaed  solutions  with  oscillating 
asymptotics  with  |p|  — »  oo.  The  spectrum  of  eigenvalues  A  of  this  set  is  continuous. 
Wth  <  1/4  it  fiUs  the  interval  (— 417*/1  —  417^,0);  otherwise  it  fills  semiaxes 
A  <  0  and  A  >  417>/(417»  - 1). 

THE  APPUCATION  OF  LUMPED  PARAMETERS  MODELS 
TO  THE  OPTIMAL  CONTROL  PROBLEMS  OF 
STRONG  BLAST  WAVES 

Optimal  contrcd  problems  of  ei^loeion-induced  shock  waves  consist  in  search  for 
optimal  kinetics  of  energy  discharge,  which  provides  the  "best”  characteristics  of  the 
shock  wave  at  given  energy  and  c^ain  restrictions  of  the  energy  discharge  rate.  The 
general  formulation  of  them  problems  is  extremely  comfdicated  both  for  blast  waves 
in  porous  dasU^lastic  media  and  gases.  Here,  we  shall  restrict  our  consideration  to 
the  most  simple  case  of  strong  explosion  in  an  ideal  gas. 
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A  significMit  nmi^fication  of  the  ptoblem  maj  be  aduered  due  to  ^>i»fication  of 
the  method  pwyoaed  bp  G.  G.  Cherapi.  The  estence  of  that  method  Hea  in  the 
uniConn]^  mfid  of  dpaamk  fonetiona  in  aeriea  over  amafl  parameta 

(t  ~  l)/(7  +  ^)>  where  7  -  is  iaentrofnc  exponent.  However,  the  firat-order  equationa 
map  be  obtained  in  aimide  qnaHtative  teaaona. 

The  modd  ia  baaed  on  the  fact  that  ahnoat  aO  gaa  maaa  behind  the  atrong  bUat 
wave  front  ia  concentrated  in  the  thin  laper  adjacent  to  the  front.The  maaa  vdodtp  V 
variea  inaignificanttp  over  the  laper  thickneaa  and  coincidence  with  the  maaa  vdodtp 
on  the  diock,  while  the  intemd  domain  preaanre  Pe  ia  levded  off  and  dependa  only 
on  time.  If  the  laper  thickneaa  to  be  neglected,  the  apatem  oi  ordinary  differential 
equationa  for  Pt  and  front  parametera  can  be  obtained  from  the  integrd  conaervation 
lawa: 


=  4ir  fl*Pe,m  = 

dt  3 


(11) 


Here  po  u  undisturbed  gas  density,  R(t)  is  front  radius  of  blast  wave,  £(t)  =  £[H(t)] 
ia  the  total  energy  value. 

FHirthermore,  we  can  apply  the  act  of  equationa  (11)  to  give  the  firat-order  equation 
for  the  value  ui  =  Jl*H*  (A  ia  shock  front  vdodty): 


3(37-1) 
7  +  1 


tti  = 


9(7^  -  1) 

4irFo 


E(R) 


(12) 


Equation  (12)  enables  one  to  determine  shock  front  parameters  and  integral  values  of 
thermal  ^Uld  kiaetic  gas  energy  for  the  arbitrary  control  function  E^R).  The  evident 
boundary  condition  for  equation  (12)  is  the  boundedness  of  ui  with  R  — »  0. 

An  extraordinary  aimpUdty  of  equation  (12)  makes  it  possible  to  reveal  the  main 
qualitative  peculiaiities  of  optimal  control  problems.  In  particular,  the  controllability 
of  blast  wave  parameters  undn  the  given  restrictions  to  kinetics  may  be  easily 
evaluated  in  that  modd.  The  possibility  of  kinetics  restoration  by  measuring  of 
shock  front  parameters  fdlows  directly  from  equation  (12). 

The  use  of  Pontrya^’s  maximum  prindple  allowa  qualitative  investigation  of  applied 
problems.  For  example,  there  has  been  performed  the  investigation  of  the  problem 
on  the  most  uniform  heating  of  finfite  gas  mass  by  the  blast  wave.  It  ia  of  interest  to 
compare  the  results  obtained  within  the  context  of  gas  dynamics  equations  with  that 
of  modd  under  discussion.  Such  comparison  was  conducted  by  the  example  of  the 
problem  of  thermal  energy  maximisation.  It  has  been  detected  that  the  results  are 
in  a  good  agreement  between  each  other,  though,  in  this  case,  the  optimal  kinetics  is 
outdde  the  formd  limits  of  the  modd  applicability. 
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CONCLUSION 


Here,  tlw  object  wea  to  demonatnte  the  poteatialitiea  and  pecuHaritiea  of  aajmptotk 
for  quaUtatire  and  quantitative  analysit  of  impact-  and  caq^lodon- 
wave  proceaaea.  The  ptobfema  diacoaaed  above  ahould  be  treated  aa  the 
exam{dea.  An  impoctant  point  ia  that  methoda  under  eonnderation  may  be  effective 
in  aolving  of  ^plied  proUema. 


BIOSKETCH 
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EXPERIMENTATION  AND  MODELING  OF  COMPOSITE  MATERIALS 


Shun-dim  Qmhi 

Materials  Krectorate 
U.S.  Aimy  Research  Laboratny 
Watotown,  MA  02172-0001 


The  discussim  will  be  focused  on  dynamic  behavior  of  continuous  fiber  reinfinced 
polymeric  cmiqwsite  materials.  Ejqterimental  techniques  using  ultrasonic 
measurements,  plate  inqract  and  high  loading  rate  pneumatic  test  machine  to 
duuacterize  conqwsite  materials  will  be  briefly  discussed. 

The  test  results  of  a  S-2  glass  reinforced  plastic  laminate  obtained  with  the 
aforementioned  testing  techniques  will  be  presented.  A  simple  mechanical 
response  material  model  will  ^  formulated  based  on  die  test  results,  and  be 
implemented  into  a  hydrodynamic  computer  program  for  simulating  a  particular 
inqtact  experiment.  Technical  issues  related  to  die  simulation  and  material 
modeling  will  be  discussed. 
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SOME  APPROACHES  TO  MODEUNGOF  COMPOSITES 
MECHANICAL  PROPERTIES  WTIH  EMPLOYMENT  OF  GRAPHS 

IHECMIY. 


V.U.Novikov 

Moscow  StHe  Open  Univenity. 

22  Korchagin  St,  Moscow,  Russia,  129805. 


1.  The  character  of  arrangement  of  composite's  structure  elements  and  their 
transference  has  big  meaning  for  medianical  prqwfties  of  matoials.  Hence  the 
I»oblem  of  structural  analysis  togedier  with  ^  problem  of  structure  modeling 
acquire  the  paramount  meaning  udmn  die  substance  syndiesis,  intqierties  and 
behaviw  under  die  external  influence  are  investigated.  Designing  die  model  by 
direct  or  indirect  mediods  it  is  possible  to  predict  die  material  structure  and  its 
numufacturing.  Perhaps  the  only  possible  and  obligatory  condition  in  this  case  is 
the  gerunetrical  equivalence  between  die  unit  model's  element  and  the  full 
material's  structure.  It  is  considered  that  die  medium  structure  units,  including 
liquid  haves  character  volume,  spatial  disposition  and  preferable  distance  between 
each  other. 

The  analysis  of  literature  on  die  composites  modeling  allows  to  select  three  main 
groups,  namely: 

•The  modeling  using  molecular  view  on  structine.  Furdier  development  of  diis 
aiqnoach  is  nano-scale  of  structure. 

•Radier  large  groiqi  of  models  with  mechano-madiematical  approaches  for 
descrqition  of  die  structure  and  properties.Some  attempts  on  using  stochastic 
apparatus  for  descriptirm  of  composite's  behavior  mechanism  was  done. 
-Phenomenological  apinoaches  have  especially  large  application.  Structure- 
phenomenological  modeling  have  die  greatest  development.  In  diis  case  die 
structure  of  conqiosite  is  described  1^  two  type  of  models: 

•type  A  -  die  set  of  components  without  consideration  of  interphase  processes; 

•type  B  •  the  same  with  type  A  ,biit  widi  interphase.  For  example,  model  A  -  two 
component  system  without  interphase,  model  B  -  three  component  system  with 
interphase. 

All  variety  of  stmcture-phenomenologica]  models  is  based  on  four  concepts  of 
cmmecting  model's  elements: 

1- successive; 

2- parallel; 

3- by  crmibmation  successive  and  parallel; 

4- using  power  function  in  the  model. 

Out  off  all  proposed  for  today  models  vdiich  ctmtain  mentioned  above  principles 
it  is  possible  to  highlight  the  following: 

A.  The  models  widi  hard  puticles  in  die  soft  matrix; 

B. The  models  with  hard  particles  in  die  rigid  matrix; 

C. The  models  with  interphase; 
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D.  Another  ones. 

te  diii  pqier  aocoidieg  to  the  graphs  theory  positions  the  analysis  of  stnne 
approaches  to  con^tostle's  modding  describe  in  litnature  for  practically  all 
known  phenomendogical  stnicnire  modek  was  done.  The  results  were  ccmqnred 
widi  oin-  own  geometrical  ^iproroh  of  modeling  the  elastic-strain  in(q)erties  of 
coa^Mdtes. 

The  piinciides  of  diis  geometrical  appixMch  has  been  presented  below. 

2.  Proceeding  from  the  siqiposttion  about  the  base  impact  of  the  fiber-matrix 
interphase  and  surface  (Hocesses  on  the  composite's  structure  and  properties  the 
agre^iqKm  model  of  conqpounded  polynm,  but  taken  as  ctunposites  bit 
("combit")  -  one-dimensional  part  of  ctmiposite  material  has  been  considered 
Mow.  Combh  is  separated  to  parts  {voportionally  wifii  ctmtent  by  volume  of 
conqmsite  con^Munds  and  the  elastic-strain  problem  is  solved. 

By  analt^  with  Om's  laws  for  circuit  the  problem  is  reduced  to  die  calculation  of 
gnqths  with  the  help  of  die  generalized  F^oph's  laws  and  the  concept  of  bar’s 
rigidity  is  introduced. 


G^=lEj/(l-M2)xSj/l.  (1) 

where  E  -  modulus  of  elasticity,  p  -  Poisson's  coeffrcient,  S  -  cross  section,  1  - 
loigdi,  i  -  combit’s  parts. 

3.  Assuming  the  hipodiesis  of  interphase  existense  in  composites  and  supporting 
the  elastic-strain  properties  continuous  from  fiber  to  matrix  using  of  the  small 
parametr's  mediod  die  formula  for  calculation  die  rigidity  of  interphase  was 
founded  as 


where  -  die  rigidity  of  matrix,  E_  and  E.  -  matrix  and  fiber  moduls,  1  - 

.  ,  ™  ,  _  .  ,  . .  ,  m  r  m 

matrix  lengdi,  o  -  mterphase  thickness. 

4.The  valiant  of  model  widi  die  normal  force  on  die  layer  conqionents  of  combit 
(x-direcdon)  was  investigated.  The  graph  of  this  model  was  drawn  and  Kirgoph's 
and  Guk's  laws  were  ajqilied  to  it 
From  (1)  and  (2)  x-mortel’s  rigidity  equals 

In  die  case  of  die  action  of  force  along  the  layer’s  components  of  combat  (y- 
directirm)  the  formula  for  y-model's  rigidity  is  eqM  to 

Gy^Gf+Cl  +  Of)^^  (4) 

where  a  =  8/31^  x  (  2  +  E^^/Ej.). 

By  die  acdmi  of  die  farce  in  Ox-direction  under  smne  angle  die  formula  is 
following: 
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(5) 


G(e  )  =  G,  G  /  V(  cos2e  +  gJ  sin2e ) 

X  y  X  y 

5.  It  is  proposed  to  consider  the  matrix  or  ciunpounded  polymer  structure  as  the 
set  of  comiiits  consisting  of  segments,  for  example  as  triai^es.  In  diis  segment 
strain  was  cormted  by  approximate  method  where  the  rigidity  G  is  described 
acconhng  to  (1). 

bi  Ae  case  of  homogeneons-polynier  matrix  medimn  it  is  siqtposed  dut  the  sides 
of  triangular  segment  have  equal  jMoperties.  For  heterogeneous  medium  die 
triangular  segment  is  designed  ^m  cmnbits  widi  parameters  of  ctHnponents  from 
dus  medium. 

Finite  element  method  was  iqiplied  for  calculation  the  triangular  segment's  strain- 
stress  state  using  polyrKunial  fimctions  for  the  description  of  defmmation. 
Rigdities  between  points  are  described  by  following  equations: 

Gj  =  E.  /(I  -  n?)  X  1/4A  X  [  +  C,  2  3  x  (1  -  / 2  (6) 

where  B  ^  .  3  and  ^^12'  gradients  matrix,  A  -  area  of  a  triangle. 

Using  gr^hs  dieory’me  corresponding  graph  is  constructed  and  its  calculation 
widi  the  help  of  Kirgoph's  generalized  laws  leads  to  receiving  the  equation  for  the 
rigidity  of  triangular  segment: 

where  digits  check  with  the  numbers  of  triangular  segments  sides  or  widi  the 
composites  components. 

In  this  piqier  ^  results  of  calculation  was  compared  widi  the  experimental 
results.  The  difference  between  diem  is  about  2-16%. 

6. The  analysis  of  die  results  of  calculation  and  specifically  die  experimental  data 
allows  to  support  that  the  "ideality"  of  the  composite's  model  is  necessary  to 
qqiroach  to  tte  "reality"  of  matorial.  The  findier  search  will  be  directed  to  the 
develo|Hnent  of  die  model  insertion-block  for  the  general  composite  model  widi 
accounting  a  disparity  between  "ideality"  and  "reality". 
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Wi  THE  ENTROPY  DAMAGE  ACCUMULATION  MODEL  OF 
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INTRODUCTION 

The  micromechanical  model  of  die  accumulated  disperse  damages  for  composites 
is  presented.  As  a  measure  of  the  damage  effect  die  increment  of  eniropic  density 
is  proposed.  The  entropic  formulation  allows  to  account  die  mechanic^  and  non- 
medumical  factors  of  fiiiliire  process. 

The  formulation  of  the  problem,  the  basic  difficulties  and  some  qiproaches  of 
modelling  and  the  bulk  of  experimental  data  are  given  in  die  papers[l-7]. 

Here  the  measure  of  local  microdefect  is  assumed  to  be  the  associated  increment 
of  die  local  entropic  density.  The  generation  and  growth  of  defects  is  modelled  in 
real  time  t.  The  variable  t  may  be  the  number  of  loading  cycles,  length  of  loading 
arc. 

The  inaease  of  entropic  density  is  siqiposed  to  be  process  accounting  the 
microdefect  eneration  and  growth,  viscous  resistance,  diffusion,  chemical 
reactions. 

THE  MODEL  FORMULATION 

The  velocity  of  entropic  density  generation  defined  by  mechanical  and  non- 
mechanical  factors  is  written  based  on  die  Coury  principle  and  Honsager  relations 
.The  first  invariant  of  die  deformation  rates,  die  gradient  of  the 
temperatore,deviator  of  the  greep  deformation  velocities  and  so  on  are 
diermodynamical  forces  of  various  nature  (scalar,  two-rank  tensor  etc.).  The 
associated  diermodynamical  flows  are  of  the  independent  factors  of  the  same 
nature. 

For  modelling  of  die  processes  it  may  be  written  the  system  of  die  government 
equations  (Honsager  relations)  that  are  connected  thermodynamical  forces  to 
flows.  Taking  into  account  the  heat  conductivity  coefficient,  concentration  of  the 
given  element  for  diffusion,  the  deviators  of  die  creep  deformations  and  stresses 
&e  increment  of  entropic  density  defined  by  viscous  resistance  and  diffusion  are 
written . 

The  one  of  die  entropic  models  may  be  formulated  widi  the  aid  of  the  energy 
balance  equation,  i^ere  member  with  of  die  fiee  energy  density  accounts  die 
entropy  increasing  based  on  the  microdefect  generation  and  growdi.  The 
distribution  of  stiffiiess  E  defined  by  microdefects  as  a  fimction  of  time  is  to  be 
known  from  experimental  data. 
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Howevo'  (he  other  model  is  ftmiiileted  whidi  pos^s  the  adveotages  of  the 
micromechanical  qiproach  md  based  <m  die  use  micrcmiechamcal  paiametm. 
The  lumdier  of  diese  parameters  may  be  defined  directly  in  experiment  and  show 
die  clear  {dqwcal  intopretatiai. 

Let  the  variable  s(t,r)  is  a  fimcti<m  of  local  microdefect  is  defined  by  vector  r. 
Based  on  local  measure  it  may  introduce  the  vector  (for  plate  cracks)  and  die 
tensor  (fnr  p(»es)  parameters  of  the  defects.  The  local  measure  of  die  damage  may 
be  formulated  as  a  symmetrical  micromechanical  two-rank  tensor  s.  It  accounts  die 
vahie  of  the  elementary  microdefect  and  orientation  of  the  defect.  The 
microparameter  p  can  introduce  as  averaging  of  die  tensor  s  over  die  accounted 
volume.  This  vsdue  is  described  also  the  irreversible  process  and  s  function  of 
increasing  of  local  entropic  density.  The  trace  of  the  tensor  p  describes  die  change 
of  the  material  volume. 

For  the  plate  microcrack  system  the  local  measure  of  microdefect  is  vector  directed 
across  the  crack  plane  and  it  let^fii  defined  die  microcrack  lengdi. 

The  increment  of  the  entropic  density  S  in  accounted  volume  is  defined  by  the 
defect  growth  and  die  generation  of  &e  new  ones.  Using  of  die  microparameter  s 
and  macroparameter  S  we  can  introduce  value  of  defects  density. 

Let  P  is  density  of  defects  (a  number  of  defects).  Then  we  can  define  die  rate  of 
the  new  microdefixts  generation.  As  a  result  die  equation  for  macroparameter  S  of 
damage  accumulation  is  written. 

The  value  s  and  the  rate  of  die  new  microdefect  generation  may  be  considered 
accordingly  as  a  thermodynamical  force  and  a  diermodynaimcal  flux.  The  value  s 
is  date  measure  of  the  dmage  effect  (the  entropy  increment)  depended  on  the 
microdefect  initiated  at  time  t. 

The  equation  for  S  nuiy  be  written  in  the  statistical  form  if  P  is  presented  as  the 
random  variable  of  die  Poisson  process.  For  example,  in  the  case  of  function  s(r,t) 
=q^[-a(t-r)].  The  equation  for  S  may  be  written  in  the  form  of  the  Langeven 
stochakic  equation  [8]  widi  the  fluctuation  member  qv. 

The  accumulation  of  darm^es  including  die  heat  conductivity,  viscous 
deformations  and  ets.  is  modelled  with  the  aid  of  kinetic  equation  for  the  rate  of 
the  entropy  density  increment  S. 

It  may  be  pointed  out  the  two  basic  micromechanical  parameters  of  the  entropy 
density  ncrement:  the  microdefect  density  P  (die  rate  of  the  microdefect  nucleation 
v(t,S)  and  die  damage  factor  due  to  the  microdefect  s(t). 

In  generally  s  is  die  tensor  object.  Therefore  P  is  die  tensor  of  the  same  rank  and 
accounts  the  defects  and  its  orientation. 

The  types  of  defects  are  depended  fi'om  the  loading  direction  with  respect  to  fires, 
the  strength  of  resin,  fibres  and  strengdi  of  bonds  between  fibres  and  resin.  To 
describe  diese  items  die  kinetic  equations  for  microparameters  s  and  P  are  written. 
Thus  we  can  write  complete  system  of  kinetic  equations  for  die  model. 

The  two  fimn  of  CDS  for  die  layered  conqmsites  are  considered  The  first  form  is 
depended  on  die  {m^agadons  of  transversal  cracks  in  the  matrix,  and  die  second 
one  is  depended  on  ^  microdamages  in  die  reinforced  elements  (fiber  crushing). 
There  are  another  forms  of  die  microdamages,  for  example,  die  fiber  breaking  is 
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iattiated  dw  cracks  in  matrix  winch  uuter  inxqngation  are  delayered  die  matrix 
along  fibers  and  dien  are  initiated  maaocracks. 

It  was  fimnd  based  <m  the  eiqieriments  that  die  polynomial  and  exponent  functions 
are  appropriate  approximations  for  the  distriburion  of  the  macrortefect  density  for 
difiBaent  stages  accumuiatimi  inocess. 

The  degradation  of  transversal  cracks  under  quasistatic  and  cyclic  loading  is 
devehqicd  wift  deacceleration  and  strengdi  stabilization.  Fm  the  sectmd  stage  of 
die  diaracteristic  dannge  states  dm  mkrodefect  density  is  increased  with 
accrieratiofL  The  transfiMtnation  to  the  second  stage  is  depemled  on  die  stress  level 
of  the  actual  damage  stteus.  To  describe  his  sta^  we  us^  mine  conqilicate  form 
of  aiqnoximation  fm  mioo^ects  density. 

The  composite  structure  is  accounted  in  the  form  of  die  stress  reanalysis  according 
the  layer^  material. 

The  analysis  of  diis  model  have  led  to  the  following  conclusions. 

1.  The  distribution  of  die  damage  status  for  the  time  is  S-  ype. 

2.  The  stage  of  die  low  growdi  of  the  damage  effect  is  sWUer  for  the  case  of  the 
microcracks  in  matrix  directed  in  fibres  and  transversal  cracks  dian  for  the  case  of 
delamination  fiom  broken  fibres. 

3.  Under  the  fatigue  loading  the  damage  accumulation  isn't 
controlled  by  the  delayering  process  . 

4.  The  heating  of  composite  is  important  factor  of  die  defect  development  process 
and  is  stimulated  essentially  the  second  stage  of  the  damage  accumulation  process. 

5.  As  a  result  of  the  diermocycle  loading  and  multisection  process  analysis  was 
found  that  the  linear  summation  hypothesis  may  be  incorrect. 

The  two-stage  approach  of  development  of  the  microdefects  is  allowed  to  model 
sadsfsctoiily  the  processes  of  die  stiffoess  degradation  and  residual  strengdi. 

The  model  of  the  stiffness  degradation  is  based  on  the  local  equilibrium  hypothesis 
and  balance  energy  equations  for  thermodynamical  potential.  According  to  the 
formulation  of  die  model  die  variation  of  the  residual  strength  with  the  damage 
status  growth  is  iqiproximated  well  (in  qualitative  form)  die  experimental  data. 
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Efficient  application  of  composites  calls  for  a  concurrent  analysis  of  materials 
selection,  design  and  manufacturability.  The  need  for  such  concurrent  engineering 
is  more  imperative  for  composites  than  for  conventional  materials  because,  for  the 
latter,  product  geometry  and  materials  selection  dictate  the  manufacturing 
processes  that  can  be  used. 

In  the  present  paper  we  describe  a  concurrent  engineering  strategy  for  composites. 
We  identify  relational  databases  that  are  needed  and  survey  some  of  the  available 
process  simulation  models  tiiat  can  be  used  to  interrogate  manufacturability  during 
design  phase.  A  unified  model  is  then  proposed  for  simulation  of  fiber  compaction 
and  resin  flow. 
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Finite  element  method  (FEM)  is  very  widely  used  in  practice  for  numerical 
calculation  of  strengdr-stress  states  and  mechanical  behavior  of  materials  and 
constructions. 

Distinctive  feature  of  the  new  program  complex  "UNIVERS"  from  traditional 
finite  element  programs,  which  have  been  usualy  written  on  procedure  languages  ( 
such  as  ALGOL,  FORTRAN,  PLl  and  etc.  ),  is  its  object-oriented  nature.  All 
algorithms,  functions  and  mediods  of  die  program  are  constracted  pursuant  to  the 
laws  of  mathematical  abstraction  (mathematical  laws  of  transformation  of  objects). 
Hierarchy  of  objects  and  laws  of  their  interaction  was  developed  too.  It  allows, 
particulaiy,  to  expediendy  use  of  vectoral  and  tensoral  nature  of  displacements, 
strains  and  stresses,  for  solution  of  the  continuous  media  and  composite  materials, 
mechanics  and  die  geomechanics  problems. 

Programs  constructed  on  the  base  of  traditional  finite  elements  method  are 
characterized  by  a  lot  of  calculations,  unsufficient  reliability  of  a  programm  code 
at  large  volumes  the  initial  text.  Such  kind  of  programs,  written  on  procedure 
algorithm  languages  are  very  hardly  madified  and  adapted  for  other  conditions 
(programming  environments,  computing  systems  etc.). 

Above  indicated  problems  could  be  avoided  by  using  of  the  object-oriented 
approach  in  programming.  Such  approach  allows  to  create  of  the  high  level,  large 
reliability,  easily  adapting  to  different  systems  programm  product.  Taking  into 
account  above  said,  algorithm  language  C  ++,  as  die  most  suitable  was  accepted 
for  development  of  a  prognunm  complex. 

The  newest  lexical  opportunities  of  ^e  language  C++  are  used  in  the  program  .  In 
particular,  finite  elements  tamplates,  principles  of  inheritance  in  hierarchy  of  data 
classes  vectoral  representation  for  storage  base  of  multy  dimensional  arrays  was 
employed. 

The  program  complex  is  intended  for  the  decision  of  the  stress-strain  state, 
stability  and  strength  problems  ; 

-  continuous  media  and  composite  materials; 

-  geomechanics; 

-  truss  constructions. 

The  finite  element  formidation  for  continuous  medium  is  realized  isoparametric. 
The  library  of  finite  elements  includes; 

1.  One-dimensial  isoparametric  element. 
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2.  Two-diineiisial  isopuametric  elements  widi  quantily  of  node  points  frmn  4  up  to 
8,  wdiicli  describe  following  conditions  stress-strain  state; 

-  plane-straine  state; 

-  plane-stress  state; 

-  axisymmetric  state. 

3.  Sotid  elements  with  quantity  of  node  points  from  8  up  to  21. 

For  finite  elements  continuous  media  simulatiiMi  the  number  of  mathematical 
models,  allowing  to  simulate  the  nonlinear  features  deformation  of  natural  and 
constructional  materiab  was  provided.  FoUowing  models  in  the  program  complex 
are  realized: 

1.  lsotro|HC  material. 

2.  Unisotropic  material. 

3.  Ncmlinear-elastic  material. 

4.  Incremental  model  of  defomability  type. 

5.  Model  of  elastic-plastic  current  widi  criterium  of  plasticity; 

-  Fon-Mises; 

-  Drucker-Prager. 

6.  Viscoelastic  material. 

As  a  whole  die  program  permits: 

-  elastic  and  nonlinear-elastic  accounts; 

-  to  use  in  accounts  simultaneously  some  materials  with  various  characteristics; 

-  to  use  in  accounts  simultaneously  some  materials  with  various  models; 

Using  of  the  "UNTVERS"  program  die  elastic-  plastic  calculations  by  one-axes 
compression  and  extension  and  pure  shear  for  die  cell  of  periodicity  of  filled 
fibrous  composite  with  interphase  layers  vere  made. 
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Research  in  mechanics  of  composites  dealing  with  physically  based  dieories  for 
the  deformation,  damage  initiation,  damage  growdi  and  failure  of  polymer  matrix 
composite  materials,  subjected  to  complex  stress  states  and  severe  environments 
will  be  discussed.  Recent  accomplishments  in  the  Mechanics  of  Composites 
research  thrust  at  the  Office  of  Naval  Research  will  be  described. 
Accomplishments  in  die  following  areas  will  be  summarized;  micromechanics 
based  duce-dimensional  constitutive  theories  for  laminated  composites;  local- 
global  tqiproaches  for  die  analysis  of  thick  laminates;  microbuckli^  dieories  for 
unidirectional  composites;  models  for  compression  failure;  delamination  growth 
models;  impact  damage  models;  moisture  absorption  and  moisture-induced 
damage  in  composites;  and  the  effect  of  defects  on  mechanical  response  and 
failure. 
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The  shape  memory  material  (SMM)  show  the  cosiderable  promise  for  the  creation 
new  structures  and  advanced  technology.  These  materials  utilize  in  deployable 
structures,  demountable  tight  joints,  ad^tive  sistems,  geometry  changing  surfaces, 
active  elements  with  high  power  or  (and)  high  accuracy  and  so  on.  The  most 
promissing  for  such  structures  are  die  composites  with  shape  memory  filaments  or 
plies  and  flexible  matrix  (SMC).  One  of  the  number  adventures  of  such  composites 
is  the  opportunity  of  the  structural  two  way  sh^  memory  effect  [1,2], 

The  hroad  introduction  of  SMM  and  SMC  is  restrained  because  of  lack  of  the 
constitutive  equations  such  a  materials  diat  would  can  describe  the  whole  set  of  its 
unique  mechanical  properties.  The  known  phenomenological  constitutive  equation 
for  SMM  [3,4]  is  onedimensional.  The  equations  based  on  the  sbructural-analitical 
approach  [5]  and  known  micromechanical  approaches  [6,7]  are  too  complex  for 
practical  purposes  in  solution  two-  or  threedimensional  problems. 

It  is  known  that  die  reason  of  the  unusual  behavior  of  SMM  are  the  termoelastic 
martensitic  transformations  prosiding  in  such  materials  under  changes  of 
temperature  or  (and)  stresses.  This  process  comprises  die  nucleation 
(dissapearance)  and  drowing  (degradation)  the  martensitic  crystals  in  austenitic 
matrix. 

Formerly  the  new  micromechanical  approach  was  elaborated  to  formulate  the 
constitutive  equations  for  solids  of  complex  structure  die  changing  of  state  of 
witch  is  connected  widi  nucleation  and  growth  some  microcarriers  [8,9]. 
Previously  diis  method  was  applied  in  continuous  damage  mechanics  where 
microdefects  play  die  role  of  microcariers.  The  damage  accumulation  equations 
was  derived  for  damage  tensor  in  the  case  of  metall  woiking  process  and  for 
damage  axial  vector  in  die  case  of  low-cycle  fatigue  theory  [8]. 

In  the  case  of  SMM  the  role  of  microcariers  play  the  crystals  of  martensite.  The 
models  are  formulate  for  nucleation,  disappearance  and  development  these  crystals 
during  forward  and  reverse  martensitic  transformation. 

In  die  framework  of  nucleation  model  the  contribution  B--  to  the  increment  of 
phase  deformation  due  to  diis  process  is  calculated  from  equation 

y  ^  y  0  ij 

where  5ij  is  die  unit  tensor,  a'ij  -  deviator  of  stress,  p,  parameters  of  material. 
In  die  model  of  distqiearance  of  the  martensite  crystals  during  reverse 
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tnunfonnation  tensor  B-.  consists  of  two  terms  corespcmding  to  the  reverse 
and  numotmiic  B.-  siiii|pe  memoty  effect.  The  unified  model  is  fonnulate  ror 
microcaiiers  development  just  as  fm  forward  so  for  reverse  transformation.  For 
this  model  die  increment  of  microidiase  defcMmation  because  of  develc^ment  of 
single  crystal  is  proportional  to  die  current  value  of  this  deformation. 

On  Ae  foimdation  of  these  micromodeles  rme  can  deduce  widi  aid  of 
micrmnerdianical  procedure  [8,9]  the  consdtiitive  equations  for  phase  deformation 
tensor  e|j^  in  SMM: 

d6.j*'/dq  =  B..  +  aQ8./  (1) 

where  q  is  die  volume  part  of  martensite. 


g  fPSij  +  Vij  if  d^K) 

'■*  if  dq<0 

Oq,  ^  is  the  material  parameters. 

Besides  of  equation  (1)  it  is  used  die  equations  for  elastic  deformation  and  termal 
expansion.  In  diese  equations  it  takes  into  account  the  difference  between  the 
values  of  elastic  modulus  and  termal  expansion  coefficient  for  martensite  and 
austenite  states. 

Unlike  die  known  proposals,  this  model  descripts  the  most  broad  set  of  properties 
of  SMM  namely: 

•  forward  martensitic  transformation  widi  constant  or; 

•  variable  stress,  includind  complex  stress  state; 

-  oriented  martensitic  transformation; 

•  sluqie  memory  including  reverse  and  two  way  shape  memory  effects; 

-  martensitic  inelastisity; 

•  pseudoelastisity  and  so  on. 

The  constitutive  equations  is  formulated  for  finite  phase  deformation  using 
Jaumarm  derivatives  for  Henky  tensor  of  deformation  and  Cauchy  stress  tensor. 

In  the  simplest  case  this  nonlinear  problem  for  two  -  or  diree  dimensional  body  can 
be  solved  analytically  widi  aid  of  Laplas  transformation.  For  more  complicated 
problem  will  be  create  die  finite  element  computer  program. 

This  constitutive  equation  for  SMM  was  used  for  investigation  of  the  behavior  of 
composites,  containing  filaments  or  plies  fi'om  SMM  and  flexible  matrix.  Two 
example  of  such  composites  are  considered  in  diise  work.  The  diirst  is 
onedirectional  composite  beam  widi  filaments  from  SMM.  The 
filaments  in  austenitic  phase  state  are  loaded  of  constant  tensile  stress  and  cool 
dnough  the  iterval  of  forward  martensite  transformation  temperature.  These 
filaments  elongate  ^  to  die  deformation  of  8%-10%.  Then  the  filaments  are 
embedded  into  the  flexible  matrix  such  as  silicone,  polyurediane  or  diflone.  If  the 
temperature  of  filaments  in  such  a  conqiosite  increase  (for  exanqile  due  to  current) 
widiin  the  interval  of  reverse  transformation  diey  shorten  because  of  shape 
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memoiy  cfiect  The  len^  of  beam  deoease.  In  this  state  matrix  have  s(Mne 
conqnesrive  stress,  fit— have  seme  tensile  <»e.  Then  if  die  tempoature  of 
filaments  decrease,  the  filaments  el<ngate  because  of  finward  martensitic 
transfinmatkm  in  the  presence  of  tensile  stress.  By  this  means  structural  two  way 
shape  memoiy  effect  is  relevant  fm  such  a  cmiqxisite.  But  die  tmsile  stress  in 
filmients  have  to  be  suffisiod  to  restitute  the  original  size  of  conqiosite.  The 
cimdition  of  closed  two  way  shape  memoiy  effect  is  derived  for  such  conqwsite 
beam: 


z^e“0(‘-l/z) 

where 

z  =  [aQ/CQE^(l-K[l-lil/tiXEj/Ejj^)) 

Ep  E  is  the  elastic  modulus  of  filament  and  matrix,  p  -  volume  finedon  of 
fiuunents.  As  well  are  found  die  relations  between  useful  load,  temperature  and 
displasement,  stresses  and  strains  in  such  a  composit.  Similarly  was  considered 
composites  beam  with  two  outer  diin  shape  memoiy  plies  and  flexible  thick  middle 
layer.  This  beam  is  able  to  twist  if  the  temperature  in  outer  plies  changes. 
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Elastically  tailined  composite  designs  are  being  used  to  achieve  favorable 
defiormation  inodes  under  a  given  loading  environment.  Coupling  between 
deformation  modes  such  as  extensi(»-twist  or  bending-twist  is  created  by  an 
qqiroiHiate  selection  of  fiber  orientation,  stacking  sequence  and  blend  of 
materials. 

Two  main  challenges  are  associated  with  the  inq>lementation  of  elastically  tailored 
cmiqwsites.  The  first,  is  concerned  with  the  development  of  a  reliable  analytical 
model  ndiich  accounts  for  dm  influence  of  the  material's  anisotropy.  The  second,  is 
die  design  of  test  mediods  vriiich  can  accurately  measure  the  various  coupling 
associated  widi  a  given  loading. 

From  the  analytical  standpoint,  the  fimdamental  mechanism  producing  elastic 
fiulming  in  advanced  cmnposites  is  a  result  of  dieir  anisotropy.  A  review  of  die 
analytical  methods  developed  for  elastically  tailored  composites  is  presented.  A 
theory  based  on  a  variation  asynqitotical  analysis  of  2D  shell  dieoiy  is  proposed 
in  order  to  assess  die  accuracy  of  predictimis  from  previously  developed  theories. 

A  number  of  methods  for  testing  extension-twist  coupling  in  elastically  tailored 
conqiosites  are  discussed.  The  benefits  and  limitations  of  each  of  die  mediods  are 
hi^ilighted.  Quantitative  results  for  each  of  die  testing  mediods  using  a  standard 
qiecimen  are  shown  and  compared  to  determine  die  qiplicability  and  suitability  of 
each  of  the  methods. 

A  comparison  of  results  widi  previously  developed  theories.  Finite  Element 
simulations  and  test  data  will  be  presented.  Two  demonstration  tests  featuring 
bending-twist  and  extension-twist  coupling  will  be  conducted  in  order  to  illustrate 
die  physical  mechanisms  which  induce  coiqiling. 
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